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Recent work has shown that dispersion relations with Planck scale Lorentz violation can produce 
observable effects at energies many orders of magnitude below the Planck energy Af . This opens 
a window on physics that may reveal quantum gravity phenomena. It has already constrained the 
possibility of Planck scale Lorentz violation, which is suggested by some approaches to quantum 
gravity. In this work we carry out a systematic analysis of reaction thresholds, allowing unequal 
deformation parameters for different particle dispersion relations. The thresholds are found to 
have some unusual properties compared with standard ones, such as asymmetric momenta for pair 
creation and upper thresholds. The results are used together with high energy observational data 
to determine combined constraints. We focus on the case of photons and electrons, using vacuum 
Cerenkov, photon decay, and photon annihilation processes to determine order unity constraints 
on the parameters controlling 0{E/M) Lorentz violation. Interesting constraints for protons (with 
photons or pions) are obtained even at 0((E/M)^), using the absence of vacuum Cerenkov and the 
observed GZK cutoff for ultra high energy cosmic rays. A strong Cerenkov limit using atmospheric 
PeV neutrinos is possible for 0{E/M) deformations provided the rate is high enough. If detected, 
ultra high energy cosmological neutrinos might yield limits at or even beyond 0{{E/M)^). 
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The principle of relativity of motion goes all the way back to Galileo [ij , who noted that observers below decks in 
a large ship gliding across a calm sea have no way of determining whether they are in motion or at rest. Einstein's 
special relativity, which is founded on this principle, has been spectacularly successful in accounting for phenomena 
involving boost factors as high as 10^^. Moreover, the Lorentz group has a beautiful mathematical structure, and this 
symmetry powerfully constrains theories in a way that has been very useful in discovering new laws of physics. It is 
natural to assume under these circumstances that Lorentz invariance is a symmetry of nature up to arbitrary boosts. 
Nevertheless, there are several good reasons to question exact Lorentz symmetry. From a logical point of view, the 
most compelling reason is that while 10^^ is a large number, it is nowhere near infinity. There is, and will always be, 
an infinite volume of the Lorentz group that is experimentally untested since, unlike the rotation group, the Lorentz 
group is non-compact. Why should we assume that exact Lorentz invariance holds when this hypothesis cannot even 
in principle be tested? 

While the non-compactness reason for questioning Lorentz symmetry is perhaps logically compelling, it is by itself 
not very encouraging. However, there are also several reasons to suspect that there will be a failure of Lorentz 
symmetry at some energy or boosts. One reason is the ultraviolet divergences of quantum field theory, which are a 
direct consequence of the assumption that the spectrum of field degrees of freedom is boost invariant. Another reason 
comes from quantum gravity. Profound difficulties associated with the "problem of time" in quantum gravity 0, Q 
have suggested that an underlying preferred time may be necessary to make sense of this physics. Also tentative 
results in string theory Q , quantum geometry |^ , and non-commutative geometry 0, H approaches to quantum 
gravity have suggested that Lorentz symmetry may be broken in the ground state. 

Finally, there have been recent hints from high energy astroparticle physics that we may already be seeing the 
effects of Lorentz violation (although as discussed below the most recent analyses make this seem unlikely.) One 
comes from the photo-production of electron-positron pairs when cosmic gamma rays collide with photons of the 
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infrared background. Below 10 TeV the (indirectly) observed absorption of such gamma rays by this process offers 
support for boost invariance up to the boost that relates the cosmic rest frame to the center of mass frame of the 
colliding photons. (For a 10 TeV gamma ray colliding head on with a 25 meV infrared photon this yields a boost of 
10^.) However, according to some (but not all) models of the infrared background, there appears to be less absorption 
than expected for gamma rays above 10 TeV coming from the blazar Mkn 501 (located at about 157 Mpc from us). 
If true this could be explained by an upward threshold shift due to a Planck scale suppressed Lorentz violating term 
in the dispersion relation for the gamma rays . 

The other hint comes from the cosmic ray events beyond the GZK cutoff 0, on high energy protons. Ultra 
high energy protons undergo inelastic collisions with CMBR photons leading to the production of pions (the boost 
to the center of mass frame yields the figure of 10^^ mentioned above). As a result, protons above ~ 5 x 10^^ eV 
are not able to reach us from distances above a few Mpc 12]. In spite of this prediction, cosmic rays with energy 
beyond 10^° eV have apparently been observed by the AGASA experiment [13 (see also |lj| for a review on this 
issue). The nature and origin of these ultra high energy cosmic rays is unknown and several explanations have been 
proposed (see 0,^3 for an extensive review). One proposal is that Lorentz violating terms in the dispersion relation 
for the proton produce an upward shift of the threshold for pion production, allowing these high energy protons to 
reach us 113, lH llE, JflJ • Interestingly it was argued that a universal Lorentz violating deformation of the particle 
dispersion relations would be capable of explaining both the TeV gamma ray absorption anomaly and the trans-GZK 
events |20) . 

The evidence for the TeV gamma ray and GZK anomalies is not convincing at this stage, however. Indeed it has been 
argued in I2H for the former and in |2^ |2^ for the latter that the data are consistent with Lorentz invariance. 
For us therefore the most important point is just that it is possible at all that Planck scale violations of Lorentz 
symmetry could be observed or constrained by current and upcoming observations. The focus of the prese nt pap er is 
almost entirely on the constraints that can be imposed. Our work extends prior results jjlTl Ha . Il9l l20l |24L l25l \26L l27| 
in several ways: (i) combining constraints to limit parameter space of a priori independent parameters, (ii) discovery 
and characterization of the asymmetric threshold effect, (iii) characterization of upper threshold effects, (iv) extending 
analysis for threshold effects to higher order nonlinearities. A brief report on some of our results has already been 
given in (2^ . Some of these results have been confirmed in . 

In the next section we discuss our theoretical framework and list the reactions we are going to consider. In Section HTTI 
we study the kinematics of some photon-electron processes in order to determine how Lorentz violating dispersion 
affects thresholds. The details of the photon annihilation threshold analysis are worked out in the Appendix. These 
results are then used to deduce observational constraints on the electron and photon deformation parameters. Taken 
jointly these constraints severely restrict the parameter plane. Section llVl is devoted to the discussion of other possible 
interactions including hadrons or neutrinos, and in section we discuss the special case of common Lorentz violating 
parameters for all the particles. Finally we present some conclusions and perspectives in section IVll 

Throughout this paper we adopt the following notational conventions: p4 denotes a four- momentum 734 = (w, p), and 
p is the magnitude of the three-vector p. The metric signature is (+, — , — , — ). We use the energy scale M = 10^^ GeV 
to form dimensionless Lorentz- violating parameters, since it is close to the Planck energy Mp = {hcycy^^ ~ 1.22-1019 
which we are presuming sets the scale for violation of Lorentz invariance induced by quantum gravity. We often employ 
units in which AI — 1. 



II. THEORETICAL FRAMEWORK AND PROCESSES CONSIDERED 



Various approaches to quantum gravity have suggested that violations of local Poincare symmetry might occur, but 
no reliable prediction is currently availble. These suggestions range from the breaking of just the boost symmetry 
to breaking of the full local Poincare group. In this paper we study the former case since it is the minimal one for 
which consequences of boost symmetry violation can be explored. Thus we shall assume that rotation and spacetime 
translation symmetries are preserved, so that in particular energy and momentum are conserved. ^ 

Dispersion relations determine how particles propagate and, via energy-momentum conservation, how their in- 
teractions are kinematically constrained. Hence Lorentz violating dispersion relations provide a relatively theory- 
independent window into the possibility of Lorentz violating physics. In this work we explore the observational 
consequences of such deformed dispersion relations in flat spacetime, i.e neglecting gravitational effects. The conse- 
quences of such dispersion relations have also been extensively investigated in the context of the Hawking effect (see 
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e.g. [33 and references therein) and the primordial spectrum of density fluctuations in cosmology (see e.g. |34ll35l l36| 
and references therein). 

In this section we discuss our framework for parametrizing such Lorentz violating physics, as well as the processes 
through which one might hope to place constraints or to observe Lorentz violation. 



A. Theoretical framework 



A dispersion relation that is not boost invariant can hold in only one frame. We assume this frame coincides 
with that of the cosmic microwave background. As mentioned above, we further assume that rotation invariance is 
preserved in this preferred frame. Thus the dispersion relation takes the form E = E{p), where p is the magnitude of 
p. In the Lorenz invariant case we have E'^ — +p^. Effective field theory suggests that it should suffice to consider 
generalizations of this form involving only integer powers of momentum, 

00 

£;2 ^m2+/ + ^a„p". (1) 

We presume that any Lorentz violation is associated with quantum gravity and suppressed by at least one inverse 
power of the Planck scale M. For n > 3 it is therefore natural to factor out the appropriate power of M and write 
o-n = Vn/M'^~'^ where rin is a dimensionless constant that might be expected to be of order unity if indeed quantum 
gravity does violate Lorentz symmetry. For n < 3 there must in addition be another mass scale, ^, which might be a 
particle physics mass scale, in terms of which the coefficents ai.2 can be written as Oi = ai fj? /M and 02 — a2 ^Ji-IM , 
where again ai^2 might be expected to be of order unity.^ In a situation such as this, the important terms at low 
energies p ^ /i would be from and p^ . At high energies p 3> /i, the p^ term if present would dominate. If this term 
is absent then the term would dominate when p^ 3> jiM . 

A large amount of both theoretical and experimental work has been done on the case n < 2. The most general 
framework is the "standard model extension" which includes not just rotation invariant effects but all possible 
renormalizable Lorentz and CPT violating terms that can be added to the standard model Lagrangian preserving 
the field content and gauge symmetries. Low energy observations js^ l38L l3fl| | have placed stringent limits on the 
magnitude of such Lorentz and CPT violating terms. For example in |39l | a very strong constraint of order 10~^^ from 
spectropolarimetry is provided for the electromagnetic birefringence of the vacuum in the standard model extension. 
High energy astroparticle phenomena |l7ll4Clj have also been used, however in the case of such phenomena the above 
discussion suggests that unless the p'^ term is absent it would be expected to dominate over the p^ and p^ corrections. 

In this paper we focus on the constraints that can be obtained from high energy phenomena. In the absence of 
peculiar tuning of the coefficients of the terms with different powers p", it is natural to suppose that the lowest nonzero 
term with n > 3 will dominate at these energies. Hence, for simplicity, we shall include only one Lorentz violating 
power of momentum. Our study thus amounts to studying the observational consequences of dispersion relations of 
the form 

^2 =p2^^2^^^p"/M"-2. (2) 

The subscript a denotes different particles, and a priori all the dimensionless coefficients rja could be different. (For 
notational uniformity we use here ryi_2 rather than the coefficients ai^2 defined above.) We assume that, in addition to 
being conserved, energy and momentum add for composite systems in the usual way."^ It might seem that the effects 
of such deformations of the dispersion relation could be important only near the Planck energy. However, there are 
at least two types of phenomena for which this is not the case. 

First, for particles that propagate over cosmological distances, small differences in propagation speed can build up 
to detectable time-of-flight differences. Second, thresholds for particle reactions can be shifted, and thresholds can 
appear for normally forbidden processes. These threshold effects can occur at energies many orders of magnitude 
below the Planck scale. To see why, note that thresholds are determined by particle masses, hence if the p" term is 



^ Renormalization group arguments might suggest that lower powers of momentum in Eq. will be suppressed by lower powers of M . 

tlowever this need not be the case if a symmetry or other mechanism protects the lower dimension operators from Lorentz violation. 

See e.g. 1371 for an example of this in a brane-word scenario where there is Lorentz invariance on the brane but not off the brane^ 

^ Note h owever that there have been recent proposals in which the composition law for energy and momentum is also modified |4ll |42L 
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comparable to the term in |5| one can expect a significant threshold shift. This occurs at the momentum 

Pdcv^ KM^-Vr,)'/", (3) 

which gives a rough idea of the energies at which we expect to see deviations from standard physics. The typical 
scales for some different particles if 77 ~ 1 are summarized in Tabled 

TABLE I: Typical energies at which one can expect deviations from standard kinematics for different particles if 77 ~ 1 and 
n = 3,4. The mass of the neutrino is taken to be ~ 1 eV, this being the current upper bound on the mass of the lightest 
neutrino. 



n 


Pdcv for i^c 


Pdcv for e 


Pdcv for p+ 


3 


~ 1 GeV 


~ 10 TeV 


~ 1 PeV 


4 


~ 100 TeV 


~ 100 PeV 


~ 3 EeV 



B. Viability of theoretical framework 

Before considering the observational constraints, a few comments are in order regarding the viability of the theo- 
retical framework we are adopting. 

• Restriction to p <^ M and monotonicity of E{p) 

We view the dispersion relation just as the initial terms in a derivative expansion, so we are assuming nothing 
about the actual Planck scale physics. In particular, when n > 2 and 77 is negative, the right hand side of the 
dispersion relation ^ becomes negative for large enough momenta ~ |r7|~^/*^"~^^M. However, we never use 
the dispersion relation in this regime where the energy would be imaginary. Moreover, it will be important for 
our threshold analysis that we restrict attention even further to the regime in which the dispersion relation is 
strictly monotonic. As long as |ry| is not much larger than unity this will be the case provided the momentum is 
below the Planck scale. In fact, we consider only momenta many orders of magnitude below the Planck scale. 

• Causality and stability 

For positive i] the propagation is super luminal at high energies. One might worry that this would lead to causal 
paradoxes, however this is not the case, since the propagation is always forward in time relative to the preferred 
frame in which the dispersion relation is specified. For negative rj the 4-momentum is spacelike at high energy, 
hence in a boosted frame the energy can be less than zero. One might think this implies that the case with 
ry < is not energetically stable and hence unviable. This is not so, however, since all energies remain positive 
relative to the preferred frame, which is enough to guarantee stability.^ 

• Dispersion relations for macroscopic systems 

The deformed dispersion relations are introduced for elementary particles only; those for macroscopic objects 
are then inferred by addition. For example, if N particles with momentum p and mass m are combined, 
the total energy, momentum and mass are i?tot = NE{p), Ptot — Np, and Mtot — Nm, so that iS^ot = 
MjQt + PjQt + N'^~"r]P^^ (in units with M = 1). The ratio of the Lorentz violating term to the term is the 
same as it is for the individual particles, r]p"~^, hence there is no observational conflict with standard dispersion 
relations for macroscopic objects. 

• Effective field theory and compatibility with general relativity 

There is no difficulty exporting deformed dispersion relations to curved spacetime, provided they can be produced 
by an effective Lagrangian for a field. In this case, the preferred frame in which the dispersion relation holds is 
specified by a unit timelike vector field, which must be promoted to a dynamical field of the theory if general 
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covariance is to be preserved |46l 1471 l48l |49( . In the cases that n is even, there are obvious Lagrangians that 
produce the dispersion relation. For example one can add terms involving extra powers of the spatial Laplacian, 
such as (*^'^)V^0)^ for a scalar field. For odd n there seems to be no local action that will work for real scalar 
fields, although for a complex scalar the term i<j)dt^^^V^(f> + h.c. induces cubic and higher order terms. To induce 
cubic terms for spinors one can write for example ?/j('^)V^?A, and for the electromagnetic field one can write 
B • V X E (which violates parity). This last case yields a sort of Lorentz violation that emerges from quantum 
geometry calculations 0. The Lorentz violating terms in the effective Lagrangians just discussed have mass 
dimension greater than four so are not renormalizable. This is not a fundamental problem, since we only regard 
the Lagrangian as an effective one below some large energy scale, however it raises the question of naturalness. 
For now we take the point of view that there may be an explanation for the low energy Lorentz symmetry that 
is not yet understood. 

C. Processes considered 

In order to determine the strongest joint constraints on the a priori independent coefficients rja in ^ one must 
identify several processes involving the same types of particles. We focus most of our attention on the case of photons 
and electrons, since the electron mass is light and these particles interact readily. In this way we are able to obtain 
rather strong constraints on the allowed parameter space. We also consider several other processes some of which 
presently allow or will soon allow further interesting constraints to be placed. Here we summarize all the processes 
to be considered in the paper and a few more. 

A ) Photon electron processes 

(a) QED vertex interactions: The basic QED vertex involves one photon and two electron lines. With all 
particles on-shell this vertex is forbidden (for any in-state) by energy momentum conservation in the usual 
Lorentz invariant theory, but it can be allowed by Lorentz violating dispersion. In particular we consider 
the following processes: 

i. e~ — > e~ 7: This vacuum Cerenkov effect is extremely efficient, leading to an energy loss rate that goes 
like E"^ well above threshold. Thus any electron known to propagate must lie below the threshold. We 
shall also discuss the vacuum Cerenkov effect for other charged particles and even for neutral particles. 

ii. 7 — > e"*" e~: The photon decay rate goes like E above threshold, so any gamma ray which propagates 
over macroscopic distances must have energy below the threshold. 

iii. e+ — > 7 Pair annihilation to a single photon can also occur. For cosmological observations this would 
be hardly distinguishable from the similar two-photon pair annihilation and as such it is not presently 
helpful in providing observational constraints. 

(b) 77 —* e+e^: Photon annihilation occurs in ordinary QED above a certain threshold, however Lorentz 
violating dispersion can modify this threshold in observationally interesting ways and can introduce an 
upper threshold. (The related reactions of pair annihilation (into two photons) and Compton scattering are 
also modified, however these effects offer no clear signal that can provide useful constraints.) 

(c) 7 ~^ Nj: Photon splitting is allowed by energy momentum conservation in Lorentz invariant QED if all 
N + 1 photons have parallel momenta, but the process does not occur both because the matrix element 
vanishes and the phase space volume vanishes. With modified dispersion the photon four-momenta are no 
longer null (and there may be additional Lorentz violating operators that mediate the process) hence this 
reaction can occur with a finite rate. However, we shall see that the rate is too small to be observable. 

(d) Time of flight constraints: Non-linearity in the modified dispersion relation leads to different times of arrival 
for photons of different wavelength emitted from the same event. Such differences can provide an upper bound 
on the parameter governing the amount of Lorentz violation for photons, independently of the parameters 
for other particles. 

(e) Vacuum birefringence constraints. Violations of Lorentz invariance involving also parity violation can lead 
to unequal speeds of propagation for different photon polarizations. The absence of such birefringence effects 
for light (IR-UV) from cosmological sources has been used to provide constraints of order 10^'^^ and 10"^ 
for the quadratic 39] and cubic deformations |5Clj| respectively. 

B ) Other processes: 

(a) Alternative vacuum Cerenkov effects 
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• p"*" — !■ p"*" 7 or n — > 717: Note that to properly analyze this reaction the structure of the proton or 
neutron must be taken into account. 

• J/ — !■ 1/7: Although neutrinos are neutral, they still have a charge structure in the standard model so can 
in principle produce vacuum Cerenkov radiation via the charge radius coupling. Massive neutrinos could 
also radiate via the magnetic moment coupling [sij . The related process of photon decay to neutrinos, 
7 — > Du^ may also provide an interesting constraint. 

• Gravitational Cerenkov radiation will occur if matter moves faster than the phase velocity of gravitons 
in vacuum (52]- This effect has been used, in the special case n = 2, to place limits on the difference 
between the maximum speeds of propagation for gravitons and photons |5,1| . 

(b) p'^ 7cMB ^ T^o- GZK interaction. Lorentz violations can change the allowed range of energies for this 
reaction. The confirmation of the standard GZK cutoff can therefore provide interesting constraints even in 
the case n — A due to the tremendously high energy of the most energetic cosmic rays. Moreover the highest 
energy events recorded by AGASA may conceivably be explained via an upper threshold. 

(c) Neutron stability -Proton instability. If the dispersion relations for the neutron and proton are independently 
modified, it is possible to make neutrons stable at high energies. The highest energy AGASA events could 
be understood in this manner if the trans-GZK particles were actually neutrons, hence suppressing their 
interaction with the cosmic background radiation. 

(d) Neutrino oscillations. Non flavor-diagonal Lorentz violations can produce neutrino oscillations, even for 
massless neutrinos j54j . For quadratic deviations in the dispersion relation (n = 2) the constraints from 
current observation have been considered in [itI Is3 . |^ leading to a constraint on the difference of speed 
between electron and muon neutrinos of about 10~^^. Constraints for higher order Lorentz violations have 
been discussed in 

(e) Anisotropy effects. The motion of the laboratory with respect to the preferred frame can produce anisotropic 
effects. Limits for the case n = 2 are discussed in [s^. Such an effect has recently been used [s^l to show that 
the Lorentz violation suggested by quantum geometry calculations is in conflict with current observations 
in Hughes-Drever type experiments. 



To obtain constraints from these reactions we shall consider the following observations: 

1. Electrons of energy up to ~ 100 TeV are inferred via X-ray synchrotron radiation coming from supernova 
remnants [58ll59l|. 

2. Gamma rays up to ^ 50 TeV arrive on earth from the Crab nebula [sif. 

3. Cosmic gamma rays are absorbed in a manner consistent with photon annihilation off the IR background with 
the standard threshold [l6l I2H . This inference depends on incomplete knowledge of the IR background and on 
assumed properties of the source spectrum however, so the consistency provides only an imprecise constraint at 
present. 

4. Different photons emitted by the same gamma ray burst all arrive at earth within a narrow time interval. 

5. The GZK cut off on UHE cosmic ray protons at ~ 5 • 10^^ eV has been observed |2^ (although events at 



higher energy may have been detected 

Using the electron-photon processes we find strong constraints on the allowable range for the photon and electron 
parameters for cubic order {n = 3) Lorentz violation, while the quartic case (n — 4) is only weakly constrained. Using 
UHE cosmic ray protons we obtain strong constraints even for n = 4. 



In this section we determine the thresholds for some elementary processes involving just photons and electrons. The 
fact that these always involve the same particles will allow us to combine the constraints provided by the available 
observations and to severely restrict the space of the Lorentz violating parameters. From here on we adopt units 
with M = 1, but occasionally display the M dependence explicitly. 



D. Observations 




III. PHOTON-ELECTRON PROCESSES 
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A. Kinematics of the basic QED vertex 



The processes e~ —^ e~7 and 7 — + e+e^ correspond to the basic QED vertex, but are normally forbidden by 
energy-momentum conservation together with the standard dispersion relations. When the latter are modified, these 
processes can be allowed. 

For photons and electrons the assumed dispersion relations are: 

a;2(fc) = e + ^k", (4) 
E^{p) = m^+p^ + r]p", (5) 

where we have introduced the notation S, — f]^ and rj = rje- Let us denote the photon 4-momentum by = (a;fe,k), 
and the electron and positron 4-momenta by p4 — (Ep, p) and 54 = {Eg, q). For the two reactions energy- momentum 
conservation then implies p^ — k4 + 94 and k^ = p^ + respectively. In both cases, we have 

{P4 - k4)'^ = ql, (6) 

where the superscript "2" indicates the Minkowski squared norm. Using the Lorentz dispersion relations Q and (O 
this becomes 



2 {EpUJk — pk cos 9) , 



(7) 



where 6 is the angle between p and k. In the standard case the coefficients ^ and 77 are zero and the r.h.s. of Eq. {T)) 
is always positive, hence there is no solution. It is clear that non-zero ^ and 77 can change this conclusion and allow 
these processes. 

We define a lower threshold as the minimum energy required for the incoming particle for the reaction to occur. (If 
the initial state is a two particle state, then a threshold is defined relative to a fixed energy for the "target" particle. 
Conversely, an upper threshold is defined as the maximum energy (if any) allowed for the incoming particle for the 
reaction to occur. Our analysis is based on properties of thresholds summarized in the following threshold theorem: 

Threshold theorem: If Ep is a strictly monotonically increasing function of p for p > for all particles, 
then all thresholds for processes with two particle final states occur when the final momenta are parallel. 
For processes with two particle initial states the initial momenta at threshold are anti-parallel. 

A detailed proof can be found in . According to the theorem, = at a threshold. This point has been assumed 
in previous work but was not shown explicitly and in fact is not true if Ep is not monotonic. 

Fixing 9 to be zero, all three spatial momenta are parallel, hence momentum conservation implies 5 = |q| = 
I ± (p — k)| = |p — fc|. In this case the relation {T)) becomes 



^fc" + r]p'' - r]\p - fc|" = 2pk 



Ep uj 
p k 



- 1 



In the situations of interest to us, the momentum p is relativistic, and the Lorentz violating terms are small: 



m/p <^ 1 



i{k/MY 
V{P/M) 



n-2 



< 1 

< 1 



(8) 



(9) 
(10) 

(11) 



Using these approximations and expanding the two energies in powers of the small quantities ((m/p)^ -1- rjp'^^ ^•') and 
^^(n-2) obtain 



Ep LO 

p k 



1 / m 

p2 



VP 



(n-2) 



1 + i^fc(«-2) 
2 



1 (?*.-.)^ 



(12) 



There is a subtlety about the truncation of this double expansion. If the ratio of the two expansion parameters is 
very large, it is possible that the second order term in one quantity is comparable to (or larger than) the first order 
term in the other quantity. In such cases, spurious results can be obtained by truncating both expansions at the 
same order. We shall proceed with the first order truncation of both expansions. One can check a posteriori whether 
the truncation is consistent. It turns out that this truncation is adequate for our practical purposes. In particular, 
although at very high energies our approximate threshold results will fail to be accurate, those energies are sufficiently 
high so as to be observationally irrelevant. 
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Another important point is that Eq. © originated from ^ together with conservation of three-momentum and 
hence it is equivalent to energy conservation E{p) — uj{k) — ±E{q). For the Cerenkov and photon decay processes 
— !■ e^7 and 7 — > e'^e~ we want only the upper and lower signs respectively, since the energy of all the particles 
should be positive. It will be unnecessary to impose this choice explicitly however, since the negative energy solutions 
are excluded by the approximations to be employed, as can be checked by just imposing energy conservation directly 
and using the same approximations. We indicate below how the approximations can exclude the negative energy 
solutions. 

Consider for example the vacuum Cerenkov process. Then energy conservation with a negative energy final electron 
reads E{p) = io{k) — E{q) with w(fc), E{p) and E{q) all positive. The smallest E{q) can be (within the monotonic 
regime) is to, so we must have ijj{k) > E{p)+m. Expanding, this becomes k + ^k^'^~^'> /2 > p + m^ /2p + rip^'^~^'> /2 + m. 
On the other hand, momentum conservation (in the threshold configuration) requires that k < p. This inequality 
implies that ^A:("~^^/2 > 77p("~^^/2 + to, which requires that either ^fc*^""^^ > 0{m/p), or |77p("~^^| ~ m/p, or both. 
In either case, we see that neglected terms such as (^fc^""^))^ are not negligible compared to the term (m/p)^ that 
has been kept. 

Truncating (|12|l at first order, and inserting the result in Eq. we obtain 



efc" + - V\P - = 2pk + |fc("-2) + 

Introducing the variable x — k/p, takes the form 



pa 



l-x-\l-x\ 



(13) 



(14) 



At threshold for either Cerenkov or photon decay p and x must satisfy this kinematic relation. Note that while 
we have assumed that p is relativistic, no such assumption is needed for the other two momenta q and k. This is 
important since we shall use (|14|l in cases where the momentum distribution is highly asymmetric. 



B. Vacuum Cerenkov effect: e 

The spontaneous emission of photons by a charged particle in vacuum is forbidden in Lorentz invariant physics since 
the sum of a timelike and null 4-momentum vectors cannot lie on the same mass shell as the timelike 4-momentum. 
Modifications of the dispersion relations of the form Q and (jSJ can allow some phase space for this reaction to 
happen. If the reaction is allowed the rate of energy loss for the case n > 2 well above threshold is dE/dt ~ ai?^, 
where E is the energy of the initial charged particle and a is the fine structure constant.^ The decay distance is thus 
only of order the microscopic distance 100/i?, hence the lower threshold of the vacuum Cerenkov effect must be above 
the maximal observed energy of any charged particle known to propagate. 

The lower threshold is the lowest value of the incoming electron momentum p for which the kinematic equation 
(|14|l has a solution. At threshold x = k/p must lie between and 1, since if k were greater than p the final electron 
momentum would have to be anti-parallel to the photon momentum, which is excluded by the threshold theorem 
(c/. Sect, nil A|l . The threshold therefore occurs at the value of x in this range for which the right hand side has a 
maximum. We substitute this x in Eq (|14|l to obtain the lower threshold momentum for the electron as a function of 
^, 77 and TO. That there is no upper threshold in this case is immediately obvious since the right hand side vanishes as 
X approaches one, allowing solutions with arbitrarily large momentum. 

The analysis is somewhat simplified by rewriting Eq. H14|l in terms of the new variable w = 1 — a;, in terms of which 
it takes the form 

o 

777 

— = ~^7i,(l - 7«)"-2 + ,y(w + . . . + U>"-1). (15) 

The relevant range of w is to 1. In the threshold configuration we have p = q + k, hence w = q/p. The general 
analysis of the threshold relations must be done on a case by case basis for different values of 71, however it is easy to 
derive partial results valid for any n. 



^ For the special case n = 2 the rate of energy loss is further suppressed by the difference in speeds, dE/dt ~ (c^ — c^)aE'^ , see e.g. llTl . 
In this case the decay distance depends on how close the two speeds are, which must be taken into account in deducing observational 
constraints on the parameters. 
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First consider the case where ^ is positive. Then the first term in (|15|l is negative, so if rj is negative there is no 
solution. If 77 is positive, the maximum of the right hand side clearly occurs at w = 1, where it is equal to (n — 1)77. 
Hence the threshold for the case when ^ and 77 are both positive and n > 2 is 

r 2 1 i/n 

Pth = 7 TT- ■ 16) 

[{n- 1)t] 

Since w — 1, this threshold corresponds to the emission of a zero energy photon. This is why the value of ^ is 
irrelevant, and the Cerenkov process takes place as long as 77 is positive. Indeed also for negative values of ^ the 
process takes place as long as 77 is positive (and even for some negative values — see below), however the threshold 
configuration may occur with the emission of a hard photon. 

One more general result can be established, namely that there is no threshold if 7; < ^ < 0. To see this observe that 
for w between and 1 we have w + ■ ■ ■ + tx;""^ > w(l — w)"^^, since the derivative of the Ihs is greater than unity 
and the derivative of the rhs is less than unity. Thus if 7; < ^ < the rhs of H15|) is nowhere positive in this range of 
w. In particular, there is no threshold in the case of equal negative parameters ^ = 7; < 0. 

The remaining parameter space for which we need to determine the threshold is the region ^ < and t] > ^. 



1. Vacuum Cerenkov thresholds for n = 2, 3, 4 



In the case n — 2 equation ({T3|) becomes 



m 



(17) 



hence the threshold occurs at it; 
given by 



1, with the emission of a zero energy photon, and the threshold momentimr is 



Pth 



1/2 



(18) 



It is clear from the above expression that no threshold exists for the special case ^ = 7/. This can also be seen 
directly from the fact that quadratic modifications in the dispersion relations are equivalent to constant (momentum 
independent) shifts in the speed of propagation. In the case of equal coefficients the electron and photon dispersion 
relations share the same Lorentz symmetry only with a modified speed of light, and hence the vacuum Cerenkov effect 
(as well as photon decay) cannot take place. Nevertheless we shall see that in the higher order cases {n > 3) these 
processes are allowed for equal positive coefficients. 

In the case n — 3 equation (|15|) becomes 



iv - 0^ + iv + 0^ 

pa 



(19) 



The form of the threshold relation for p depends on the values of rj and ^. We find two different formulae, depending 
on whether the threshold occurs with emission of a low energy photon {w ^ I) — which we label as case a) below — 
or with emission of a photon with energy of order p — which is labeled case b). After a bit of calculation we find: 



Pth 



2r] 



1/3 



b) Pth 



4 77^^ (^ + rj) 



1/3 



for 77 > and ^ > —877, 



for ^ < -37y < or ^ < 7; < 



(20) 



(21) 



c) No threshold 



for 77 < and £,> rj. 



(22) 

In case b) the value oi w = q/p at the threshold is given by q/p = — (ry — £)/2{r] + £). Given a maximal en- 
ergy/momentum Pmax for which no vacuum Cerenkov effect is observed, the constraint on the parameters can be 
written as: 



a) 77 < 



2p- 



,3 ' 
max 



b) ^ > 77-2 



-2\ 



-27? 



Pmax 



(23) 
(24) 
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The case in which the correction is of quartic order is similar to the cubic one, ahhough somewhat more comphcated. 
In the case n — 4 equation H15|l becomes 

9 

— = {v-0w + {v + '20w^ + {v-0w^- (25) 
P 

With the definitions A = 77 — ^ and r = (r; + 2^)/A, (|25|l takes the form rn^/p* = X{w + tvj^ + w'^), which is what 
we used to carry out the threshold analysis. Again the form of the threshold relation for p depends on the values of 
77 and ^, and we label the cases with a) and b) as for n = 3. After some tiresome analysis we obtain the following 
expressions: 

a) = [ — j for 77 > and ^ > -(8 + 6^/2)77, (26) 

b) Pth = , for e < -(8 + 6\/2)77 < or ^ < 77 < 0, (27) 



c) No threshold for 77 < and i> r\, (28) 



where the function _F(A,t) is given by 



nA,r).AA 



> 2 



(29) 



In the case a) we have again the emission of a low energy photon (w V). In case b) the value oi w — q/p a,t the 
threshold is given by q/p — (— r — v — 3) /3. So for n = A given a maximal energy /momentum (say Pmax) for which 
no vacuum Cerenkov effect is observed, the constraint for case a) can be written as: 

9 

Tfl 

a) V<ir^- (30) 



3p: 



max 



The constraint for case b) has a cumbersome form but the corresponding line in the ^-77 plane can be found from 
equation ^n\. 

2. Observations and constraints from absence of vacuum Cerenkov effect 

We can now consider the actual constraints observations impose on ^ and 77. The previous analysis shows that the 
smallness of "rr? jp^ax determines the strength of the constraint provided by the vacuum Cerenkov effect, hence the 
strongest constraint will be obtained by considering the highest energy observed for a given particle. 

Electrons in particle accelerators are stable against the vacuum Cerenkov effect at energies up to 500 GeV, and 
in cosmic rays energies of '--^ 2 TeV have been detected |2^ |63|- Even higher energies, in the range 50 — 100 TeV, 
are necessary in order to consistently explain the peaks in the X-ray and TeV regions of the photon emission from 
supernova remnants such as SNR1006 or the Crab Nebula |2^ |5^ |^. In particular for the Supernova remnant 
SN1006 a clear identification of a synchrotron emission together with the independent estimate of the magnetic field 
strength allows one to infer that the electrons should have energies of about 100 TeV [ssL Is^.^. These electrons 
propagate over distances far longer than that required by the vacuum Cerenkov effect to decrease the electron energy 
below the threshold.^ 

For n = 2 we see from Eq. (|18|) that (t; — ^) < 3 x 10^^^ which can be compared with the limit (t; — ^) < 5 x lO"^'^ 
obtained by Coleman and Glashow ilTj using a Pmax of 500 GeV. The Cerenkov emission rate (cf footnote EJ) is fast 



After this work was completed we found l6Ci that the synchrotron emission is sensitive to Lorentz violation, and in fact one cannot be 
certain about the existence of these 100 TeV electrons for positive r\. However one can instead use the existence of 50 TeV electrons 
inferred from the detection of 50 TeV photons produced by inverse Compton scattering in the Crab nebula. This would weaken the 
constraint by just a factor of 2^ = 8. 

The competing energy loss by synchrotron radiation is irrelevant for this constraint. The rate of energy loss from a particle of energy E 
due to the vacuum Cerenkov effect goes like — e^E-^, while that from synchrotron emission goes like —e'^B^E^/m'^ (using units where 
c = ?!. = 1). For a magnetic field of about one micro Gauss (such as those involved in supernova remnants) the synchrotron emission 
rate is 40 orders of magnitude smaller than the vacuum Cerenkov rate. 
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enough for such parameters that /S.E ~ E over a distance scale of centimeters. For n > 3 the emission rate is 10^^ 
times higher. For the cases of n = 3 and n = 4 the corresponding value of rn^ /p^ax is ~ 3 x 10"'^ and 4 x 10^^ 
respectively. We therefore obtain an interesting constraint for the cubic case but not for the quartic case, assuming 
that the Lorentz violation is at the Planck scale. (We shall see in section HV A II that one could get a good constraint 
even for the n — A case by considering the 10^*^ eV cosmic ray protons, modulo some caveats that we shall discuss.) 
Figure n shows the excluded region for the parameters <^ and 77 in the n — 3 case as determined by the conditions H23|l 
and if^ . 




FIG. 1: Constraint from the absence of vacuum Cerenkov effect for n — 3. The filled region in the parameter space is the 
one not compatible with the existence of the ~ 100 TeV electrons indirectly detected via synchrotron emission from supernova 
remnants [Sg. The point where the vertical line crosses the rj axis is 77 = ?n^/(2pmax) ^ 1-5 x 10"''. 



C. Photon Decay: 7 — > e"*" e 



The spontaneous decay of a photon into a electron-positron pair, is another reaction usually forbidden by energy- 
momentum conservation. As in the case of the vacuum Cerenkov effect, modifications of the dispersion relation allow 
this reaction to occur. By the threshold theorem (c/. Sect. IlII A|l . we know that the final particles have parallel 
momenta, so that both lepton momenta are less than or equal to k. Thus x :~ k/p > 1 in Eq. H14I) . so that 
|1 — a;| = a; — 1. It is convenient to use the variable y :^ 1 / x ~ p/k, whose relevant range is zero to one. In terms of 
y, takes the form 



"fc" 



e?y (1 - y) - TO (1 - M"-'^ + (l-y) 



(31) 



The threshold corresponds to the maximum of the right hand side of Eq. (|31|l with respect to y. Note that the rhs 
is symmetric about y = 1/2, since the two leptons are kinematically interchangeable, hence it is always stationary at 
y = 1/2. However, this stationary point can be a maximum or a minimum, depending on the values of 77 and ^. If it 
is a maximum the threshold momentum is given by 



«th 



2"m2 



2(n-2)^_ ^ 



l/n 



(32) 



In the special case £, — rj, which has been mostly studied in the literature, it can be shown that the only stationary 
points of (|31|l are y — 0, 1/2, 1. Given that the right hand side of Eq. H31|l is always zero at y = 0, 1 it follows that for 
equal coefficients the threshold condition is always realized with a symmetric distribution of the final momenta. 

Contrary to relativistic intuition, and to what has been assumed in all previous calculations as far as we know, the 
threshold does not always occur with the symmetric configuration. The reason is that when 77 < 0, the lepton energy 
E{p) has negative curvature E"{p) < for sufRciently large momentum if n > 2, unlike the usual Lorentz invariant 
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case. If the threshold Ues within the negative curvature region, it cannot occur with the symmetric configuration since 
the energy of the final state at fixed momentum could be lowered by making the momentum of one particle smaller 
and one larger by an equal amount. For ?/ < ^ < 0, the threshold does occur in the negative curvature region, hence 
it is asymmetric. 

The occurrence of the asymmetric threshold might seem especially surprising if wc think, with relativistic habits, 
that at threshold the electron and positron should be created at rest in the center of mass frame. The error lies in a 
misleading application of the Lorentz transformation in the case where a definite preferred system exists. First, the 
center of mass frame may not even be accessible if the photon energy-momentum vector is spacelike (i.e. subluminal 
dispersion). Second, even if we can boost to the center of mass frame, in this frame the dispersion relation of the 
electron/positron may not have its minimum energy at zero momentum. Therefore it is not always true that the final 
particles are produced at rest in the center of mass frame. 

We now examine the cases n = 2, 3, 4 individually. 



1. Photon decay thresholds for n — 2,3,4 



In the case n — 2 Eq. ()31|l takes the form 



For ^ — < there is no threshold, while for ^ 
the threshold formula 



In the case n = 3 H31|l takes the form 



/fc2 



(C-'7)y(i-y)- 



(33) 



77 > there is a lower threshold at y = 1/2. In this case one obtains 



fcth 



2m 



(34) 



m 



(1 - y) - m (i -y) + (i - y) 



(35) 



To determine the threshold we need to find the maximal values of the rhs. The task of finding the maxima is simplified 
by introducing the new variable z = (2y — 1)^, so that y = (1 + y/z)/2, (l — y) = (1 — y/z)/2, and y{l — y) = (1 — 2;)/4. 
The relevant range of 2: is [0, 1], where z = corresponds to the symmetric configuration y = 1/2 and z — 1 corresponds 
to y = 1. 

In terms of z, becomes 



P 4^ 



^{l-z') 



(36) 



The symmetric extremum at y ~ 1/2 corresponds to z = 0, and there is one other (asymmetric) extremum at Za — ^/rj. 
One of the two extrema is a maximum and the other is a minimum. Since the second derivative with respect to z is 
77/4, the one at Za is a maximum ^ if and only if 77 < 0, and it lies between zero and one in this case if and only if 
77 < ^ < 0. Note that in the special case S, = rj the asymmetric threshold solution is removed and a threshold exists 
just for positive values of 7;. 

The value of the rhs of (|36|l at z = is (2^ — 77)/8, while at z — Za it is — (^ — 7^)^/877. We thus see that photon 
decay is allowed only above the broken line in the 77-^ plane given by f = 7//2 in the quadrant ^, 77 > and by ^ = 77 
in the quadrant ^,77 < 0. Above this line, the threshold is given by 



a) k 



•th 



8rn2 



2 \ 1/3 



b) k 



•th 



2^-77 
—Srn^ri 



for ^ > 0, 



1/3 



for 77 < ^ < 0. 



(37) 
(38) 



This does not also show that the extremum at z = is a maximum, since the relation between z and y is not smooth there. In fact, 
cP/dy^ = I62 d'^/dz^ + 8 d/dz, so at 2 = we have d^/dy^ = 8 d/dz. Using this we see that the symmetric solution is a maximum if and 
only if ^ > 0, so the asymmetric solution is the maximum if and only if g < 0. This is the same condition as < 0, since if Za = 5/*? is 
greater than zero, ^ < if and only if 77 < 0. 



13 



The detection of gamma rays with momenta up to some fcmax imphes that the parameters must he in the ^ — 77 plane 
below the line corresponding to a threshold at fcmax- This translates into the following constraints for the parameters 
^ and 77 



Am 



b) ^<r] 



1,3 



(39) 



In the case 77 = 4 Eq. H31|l can again be conveniently rewritten in terms of the variable z introduced after H35|) 
above, yielding 



^ = ^(1 
4^ 



2z-3z^). 



(40) 



The asymmetric extremum here occurs at Za = (2^ + ri)/3ri. This is again a maximum if and only if 7/ < 0, and it 
lies between zero and one in this case if and only if 7; < ^ < —rj/2. Note that again in the special case ^ = rj the 
asymmetric threshold solution is removed and a threshold exists just for positive values of 77. 

The value of the rhs of of (|40|l at z = is (4^ — 77) /16, while at z = it is —{^ — -q)^ /12r]. We thus see that photon 
decay is allowed only above the broken line in the 77-^ plane given by ^ = 77/4 in the quadrant ^, 77 > and by ^ = 77 
in the ^, 77 < 0. Above this hue, the threshold is given by 



b) fcth = 



16777,^ 

— 12777^77 



1/4 



for ^ > -77/2, 



-,1/4 



for 77 < ^ < -77/2. 



(41) 
(42) 



Again, given a maximal observed momentum for which gamma decay is not observed gives constraints on the 
parameters ^ and 77 



»7 



4777^ 
max 



b) C<^- 



12771^77 

max 



(43) 



2. Observations and constraints from absence of photon decay 

We can now consider the constraint on ^ and 77 imposed by the absence of photon decay in current observations. 
As before, the smallness of rn? /k"^^^ determines the strength of the constraint, hence the strongest constraint will be 
obtained by considering the highest energy photons observed, which are the 50 TeV gamma rays arriving on earth 
from the Crab nebula [61|. The rapid decay rate (F ~ £^ above threshold) implies that in order to propagate at all, 
let alone to reach us from the Crab nebula, these photons must have an energy below the threshold. For the 50 TeV 
photons we have m'^/kJ^ax ~ ioi4"--44_ p^j. — 2 and 77 = 3 this yields strong constraints on ^ and 77, however for 
71 > 4 this number is ~ 10^^ so the constraints are not so interesting. The case ti = 2 has already been studied 
in [f^ l40l|. Reference also uses 50 TeV, which from Eq. (|34f) yields the constraint ~ rj) < 10^^^. Here we 
consider the case 7^ = 3. 

In the case 77 = 3, we use expressions (|37|l and (|38() for the threshold momenta to impose the condition that photon 
decay be forbidden for photons below fcmax — 50 TeV. This defines a broken line in the ^-77 plane below which the 
coefficients must lie: 

a) C < I + 0.08, b) ^<Ti+ V/-O.I67;. (44) 

Constraint (a) applies for ^ > while (b) applies for ^ < 0. The excluded region in the parameter space is is shown 
in Fig. 13 

The joint constraints imposed by both vacuum Cerenkov and photon decay are shown in Fig. |31 We see that these 
two reactions are already enough for ruling out most of the parameter space. Next we shall see that by taking into 
account also the process of photon annihilation this constraint can be further improved. 
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FIG. 2: Constraint from the absence of photon decay. The filled region in the parameter space is the one excluded by the 
observation of gamma rays of energies up to ~ 50 TeV. 
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FIG. 3: The graph shows the combined observational constraint derived from the absence of vacuum Cerenkov effect and 
gamma decay. The horizontal blue shading identifies the region excluded by gamma decay, the vertical red one by Cerenkov. 
Although not visible there is a tiny region of positive ^ and r/ allowed by present observations, and there is a barely visible 
region of positive ^ and negative rj. Also the diagonal is in the interior of the allowed region. 



D. Photon annihilation: 'yJ—^e'^e 

In standard QED two photons can annihilate to form an electron-positron pair. If one of the photons has energy luq, 
the threshold for the reaction occurs in a head-on collision with the second photon having the momentum (equivalently 
energy) ki^i = m? /uiq. For /cli — 10 TeV (which will be relevant for the observational constraints) the soft photon 
threshold is approximately 25 meV, corresponding to a wavelength of 50 microns. 

In the presence of Lorentz violating dispersion relations the threshold for this process is in general altered, and 
the process can even be forbidden. Moreover, as noticed by Kluzniak in some cases there is an upper threshold 



15 



beyond which the process does not occur. ^ In this section we discuss how the thresholds depend on the Lorentz 
violating parameters. We then discuss the observational consequences and constraints that can be obtained using the 
absorption of TeV gamma rays of extragalactic origin by the intervening infrared (IR) background. 

The threshold equation for photon annihilation can be obtained by modifying our previous analysis of photon 
decay. The difference is that the initial state includes two photons rather than one. We are interested in the case 
where one of the photons has low energy (IR), hence for that photon the modification in the dispersion relation 
can be neglected. The threshold theorem (c/. Sect. IIII All tells us that the threshold configuration is a head-on 
collision. Denoting the IR photon energy by wq, the total four- momentum of the initial state thus takes the form 
/C4,in = (w(fc) + Wo, fc - ujo, 0, 0) =: (w', k', 0, 0). 

To adapt our previous calculation, we need only replace k hy k' ~ k — ujq and uj{k) by ui'{k') = uj{k' -\- loq) + ujq. 
Expanding one gets u){k' + ujq) = k' + ujq + {^/2){k' + o^q)"^^ + • • •. Since Wg ^ k, and the last term is already 
Planck-suppressed (or, if n = 2, suppressed by the small value of we can neglect ujq in that term. This yields the 
approximation Lu'{k') = k' + /2){k'Y"'~^\ where ^' is defined by 

The kinematic equation for photon annihilation is thus obtained from that for photon decay 1)31(1 by the replacements 
k ^ k' on the Ihs and ^ ^ ^' on the rhs. We can further neglect the difference between k' and k on the Ihs since 
luq ^ k, hence to a sufficiently good approximation we can use the kinematic equation 

= F{k, y) + (e + 2/ (1 - - w (1 - v) [y'""'^ + (1 - vt^'^] ■ (46) 

The variable y is defined by y = p/k, where p is one of the lepton momenta. Our analysis of the thresholds is based 
on Eq. (|^ . 

As in the case of photon decay, the thresholds occur at the symmetric value y — 1/2 only for certain ranges of the 
parameters ^ and "q. The analysis for photon annihilation is more complicated however since for n > 3 the dependence 
of Eq. 1(4611 on k and y does not separate, unlike in Eg. ((31(1 . Thus it is not simply a matter of finding the value of 
y between zero and unity for which F{k,y) is maximum. Analyzing the threshold structure is a rather lengthy and 
complicated process, so we have placed the details in an Appendix. The analysis reveals a number of unexpected 
features that thresholds can have in the presence of Lorentz violating dispersion, with intricate dependence on the 
Lorentz violating parameters. Here we summarize the results in the cases n = 2,3, and apply them to obtain further 
observational constraints. 

We obtain results valid for any value of the soft photon energy too and "electron" mass m by employing appropriately 
scaled quantities: 

fi^k/ku, ^=,7(m2("-i)/a;«), ^ - ^(™'<""'V'^o) (47) 

where ^li is the standard lower threshold m? /loq The basic threshold structure will be given in terms of these 
variables. For the case of most interest to us, n is 3 and m is the electron mass. For loq — 25 meV we then have 
^ — (wq/tti"') ^ ~ 2.3^, and similarly for 77. 

It is worth noting that while we have been thinking of as fixed and determining the corresponding high energy 
threshold, it can be viewed the other way around. The parameter (3 can also be written as luq / [m? / k) . If now k is 
considered fixed then wq is the modified soft photon threshold and rn? /k = ujia is the corresponding Lorentz invariant 
threshold. Therefore (i has also the interpretation wq/wli, that is the factor by which the soft photon threshold is 
shifted at fixed hard photon energy k. This interpretation is valid for lower thresholds only however. There is in fact 
never an upper threshold for the soft photon at fixed k (as long as loq <^ k). 

1. Photon annihilation thresholds for n=2 

For n = 2 the threshold configuration is always the symmetric one. The contour of threshold /? is given by the 
straight line 

1-0 

^ = ^+4—f. (48) 



As discussed below in section fill D 31 our results agree with those of l2i^ only in certain limiting cases. 
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The ^-intercept decreases monotonically from cx3 to — 1 for /3 < 2, and increases monotonically from — 1 to for /? > 2. 
Hence the process is forbidden below the hne ^ = ^ — 1. The parameter j3 gives the lower threshold for (3 < 2 and the 
upper threshold for /3 > 2. If the lower threshold is greater than unity, then the upper threshold exists and is given 
by P/{(3 — 1). The maximum lower threshold (3 — 2 corresponds to fc = 2kiA = 2rn? /loq. 



2. Lower threshold of photon annihilation for n — Z 



For n = 3 the threshold configuration is not always symmetric in the outgoing momenta. Instead of straight parallel 
lines for the contours of threshold [3 we find a more complicated structure. FigureQlshows the regions in the parameter 
plane where the threshold configuration is symmetric or asymmetric or does not exist at all, and a contour plot of the 
lower threshold is shown in Figure |S| 
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FIG. 4: Regions where the lower threshold for photon annihilation with n = 3 is determined by the symmetric configuration 
(light grey region), the asymmetric one (dark grey region) or the reaction does not occur (white region). The dotted fine is the 
locus of points where the contour of constant {3 < 1.5 switches smoothly from the asymmetric to the symmetric solution. 



The threshold can be symmetric only for [3 < 1.5. The symmetric part of the contour is given by the straight line 

restricted to the region above the line ^ = — 4//3^. Below this line the /3-contour switches to the asymmetric threshold, 
and is given by 




i = ^-^ + \-l^- (50) 



The joining point of the symmetric and asymmetric parts of the /3-contour is at (??join, Cjoin) = {—8/f3^, —4//?^). As 
(3 varies from to 1.5 these joining points trace out the curve ^join = ~{~'>f)'^^^- The asymmetric threshold contours 
for f3 > 1.5 terminate at the symmetric P = 1.5 contour, and accumulate above the diagonal as /3 — s- oo. The precise 
degree of asymmetry at threshold, i.e. the ratio of electron momentum to incoming hard photon momentum, is given 
by 2/ = (1 ± ^)/2, where Za = + 



3. Upper threshold of photon annihilation for n = 3 



Upper thresholds exist for n = 3 only below the diagonal and between the /? = 1.5 and j3 = oo (which gives the 
same line as (3 = 1) symmetric contours (|49|l . For a given (3 the threshold is symmetric in the region above the 
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FIG. 5: Contours of constant lower threshold /3 = const. For /9 < 1.5 these correspond to symmetric configurations above the 
dotted line and asymmetric ones below. For /3 > 1.5 there are only asymmetric lower thresholds and the contours terminate 
at the symmetric /9 = 1.5 contour. The end of the dotted line is at (??, 5) = (—64/27, —48/27), and the /? = 1.5 contour meets 
the diagonal at = ^ = —32/27. 

line ^ = — 4//3^ and asymmetric below, where the contour is given by the curve (|5()|l . The regions of symmetric and 
asymmetric upper thresholds for n = 3 are shown in Figure El The boundary of the lens shaped region next to the 




FIG. 6: Regions where the upper threshold is determined by the symmetric configuration (light grey region) or the asymmetric 
one (dark grey region). In the white region below the light grey and below the diagonal the reaction never occurs, and in the 
rest of the white region there is a lower threshold but no upper threshold. 

diagonal is determined by the curve ^join = ~(— J?)^^^ consisting of the points where the symmetric and asymmetric 
segments join. The bottom of the lens meets the diagonal at = ^ = — 1 where the symmetric ji — 2 line crosses, 
so asymmetric upper thresholds exist only for /3 > 2. The lower boundary of the region of upper thresholds is the 
/3 = 1.5 line, which meets the diagonal at 77 = ^ = —32/27. 

The possibility of upper thresholds for photon annihilation has been previously discussed by Kluzniak ||25|], who 
gave results for the values 77 = 0, ^ = — 1, and = ^ = — 1 in the n = 3 case. It seems that only the symmetric 
configuration was examined in 25], hence his results cannot fully agree with ours in cases where the asymmetric 
configuration is important. For the case 77 = 0, and negative ^, our results show that there is a symmetric upper 
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threshold only for ^ values above the P = 1.5 line, i.e. for ^ > —16/27. Our upper threshold agrees with that of 
in the limit |^/4| = |^m''/4u;Q| ^ 1. The left hand side is unity for ^ = — 1 and luq ~ 20 meV, hence our results agree 
approximately provided luq is greater than about ~ 40 meV. In the diagonal case, while our results for the symmetric 
configuration agree in the same limit, we have seen that there is no upper threshold since asymmetric configurations 
exist for arbitrarily large (3. 

4- Observations and constraints from absence of deviations from standard photon annihilation 

The Cerenkov and photon decay constraints leave open an infinite wedge-shaped region including the diagonal in 
the lower left quadrant for the case n — 3. A constraint from agreement with standard photon annihilation would 
be complementary to these and hence has the potential to confine the allowed region to a small neighborhood of 
the origin. Such a constraint is provided by indirect observations of annihilation of high energy gamma rays from 
blazars on the cosmic background radiation (CBR). Since there is presently considerable uncertainty regarding both 
the background radiation and the nature of the sources, the constraint that can be extracted is not yet very precise 
however. 

Another limitation of the present work arises from the fact that each observed gamma ray has the opportunity to 
interact with soft photons at any energy above the threshold, so to compare with observation one should compute the 
absorption using the Lorentz violating dispersion relation, integrating over all target frequencies. Such an investigation 
lies outside the scope of the present paper, so we shall only attempt to roughly characterize how large a threshold 
shift might be compatible with current observations. 

We now summarize the observational situation. The BL Lac objects Mkn 421 and Mkn 501 are a type of blazar 
emitting high energy gamma rays whose observed spectrum reaches 17 TeV in the case of Mkn 421 \m and 24 
TeV in the case of Mkn 501 '6(y\. The source power spectra are reconstructed accounting for absorption via photon 
annihilation on the intervening CBR, which ranges from the near infrared (NIR, ~ 1 fun) to the cosmic microwave 
background (CMBR, ~ 1000 /xm). Currently we have a good knowledge of the NIR and CMBR but uncertainties 
remain regarding the distribution in the intermediate, mid infrared (~ 10 /im) and far infrared (100 /im), regions 
(see e.g. Figure 1 of |63| or the discussion in llj). Some models of the IR background imply a source spectrum 
for Mkn 501 with an unexpected amount of radiation (a "pile-up") above 10 TeV 1^ such IR backgrounds 

are correct, the pile-up might be due to a process producing enhanced emission at energies larger than 10 TeV [6^ . 
or it might be explained by anomalously low absorption caused by an upward shift of the threshold due to Lorentz 
violation "20", '25, "2^, '27^ . However, recent work 21, 4^ based on improved reconstructions of the FIRB and on a 
new analysis of the gamma ray flux from Mkn 501 supports the view that current observations are consistent with the 
predictions of standard Lorentz invariant theory up to 20 TeV. Even without resolving the question of the pile-up, it 
seems well established that some degree of photon absorption has been observed up to 20 TeV, which already provides 
an interesting constraint on Lorentz violation. Moreover, it is our impression that the suggestions of an anomaly above 
10 TeV will likely prove illusory as new observations are made available, confirming the results of 0,113 We can 
thus obtain observational constraints from the requirement that the Lorentz violation does not too strongly modify 
standard Lorentz-invariant thresholds for photon annihilation. The strength of the constraints depends of course on 
the order n of the Lorentz deformation. The general threshold equation HA3|I shows that an order unity constraint 
on /3 translates into an order unity constraint on rj and ^, which corresponds to an order /m^^""^^ constraint on 77 
and ^. Since all studies seem to agree that more or less standard Lorentz-invariant absorption is occurring for gamma 
rays up to 10 TeV, we shall use the corresponding soft photon threshold of loq = 25 meV ~ 50 fj-m as a numerical 
benchmark. One then has Wg/m^ ~ 10~^^ for n = 2, Wp/m^ ~ 1 for n = 3, and Wq/ttt.^ ~ 10^^ for n = 4. Hence only 
the n = 2 and n = 3 cases can provide interesting constraints. Note that in the n = 3 case, which is of most interest 
to us, the dependence on ujq is cubic, so for example a constraint at 2 luq is eight times weaker than a constraint at 
Wo J while one at ljo/'2 is eight times stronger. This means also that there could be strong deviations in absorption 
for, say, 20 TeV gamma rays, and yet little deviation for 10 TeV gamma rays, since the standard soft target threshold 
w? /E is half as large for the 20 TeV gamma rays. 

To formulate the constraints we begin by identifying the contour in the ^-77 plane, for which the threshold is 
not shifted away from the Lorentz-invariant value. For n = 2 this no-shift contour is given by the diagonal £, = rj 
(corresponding to equal speeds of light for electrons and photons), which is independent of the soft photon energy loq. 
For n = 3 the contour is given by the joined symmetric and asymmetric 13—1 contours ()49|l and (|50|l converted to 



After this work was completed a furt her observational analysis appeared |6Sll . This allows the observational basis for the constraint 
discussed in this paper to be solidified ■ 
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the unsealed parameters, 



2 



for T] > — 8tJo/m'' 



(51) 




otherwise 



(52) 



The symmetric part is independent of luq but the joining point and the asymmetric part are not. 

Above the no-shift contour, Lorentz violation lowers the threshold. Since the shift would be larger for higher 
energy gamma rays this might, depending on the details of the IR background spectrum, enhance the "pile-up" in 
the reconstructed source spectrum if the IR backgrounds of 9] are used, or it might produce a pile-up where one did 
not otherwise exist if the IR background of jU] is used. We thus consider it unlikely that there is much downward 
shift of the threshold. In any case, nearly all of the region above the no-shift line is already excluded by the photon 
decay and Cerenkov constraints. 

Below the no-shift contour, Lorentz violation raises the threshold. We now consider the constraints this can yield 
in the cases n = 2 and n = 3. 

a. n = 2 Photon annihilation constraints. Constraints in the n = 2 case have been previously examined in Ref. 
|40|. although it was not realized there that the maximum upper shift is P — 2, beyond which the process does not 
occur at all. The /3 — 2 contour 148|) is a line of unit slope and ^-intercept —1 in the scaled parameters, hence unit 
slope and ^-intercept — cjg/m^ ^ —10^^^. As long as the 25 meV photons annihilate at least with 20 TeV photons 
(whose normal threshold is 12.5 meV), the parameters must lie above this line. 

b. n = 3 Photon annihilation constraints. For n = 3 the contours of constant threshold in the scaled parameters 
rj and ^ are shown in Fig.|Sl The process does not occur for parameters below a broken line consisting of the diagonal 
up to ^ = 77 X m'^/wg = —32/27, and the line of slope 1/2 for greater rj. If absorption at ujq is occurring for any hard 
gamma ray, the parameters must lie above this broken line, so in particular everything on and below the diagonal is 
excluded for rj < —32/27. For ojq = 25 meV this corresponds to 77 < —2.3 • 32/27 « —2.7. This is important, since 
it is a strong constraint excluding most of the diagonal, which has been preferred by some researchers |2(1 l27j . It 
is likely that a much stronger constraint holds however, restricting the lower threshold at 25 meV to be not more 
than some number of order unity times its usual value. We have indicated in Fig. |7|the form of the region below the 
no-shift contour and above the shift-less-than-/3 contour for f3 equal to 10, 5, 2 and 1.5. A stronger constraint would 
not exclude more of the diagonal, but it has the potential to chop off the infinite wedge of Figure |3| at around the 
same place it excludes the diagonal. 



FIG. 7: The unfilled region indicates parameters allowed if the lower threshold for a soft photon of 25 meV is (a) not shifted 
down and (b) not shifted up by more than 1.5, 2, 5, 10, and infinity. The upper line is the no-shift contour. No curvature due 
to the asymmetric solution is visible for this line because the junction point as defined in Eq. I5H is at 77 = —20. The line 
for the existence of a lower threshold is the lowest line. It is coincident with the symmetric (3 — 1.5 line below the diagonal 
and with the (dashed) diagonal below the crossing point. The curves stemming from the /3 = 1.5 contour are the asymmetric 
contours for /3 = 10, 5, 2, with lower values of /3 corresponding to the curves with less slope. 
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E. QED processes without thresholds 



We now consider two QED effects that occur in the presence of Lorentz violation without any threshold, velocity 
dispersion of photons in vacuo and photon splitting. The former will eventually provide competitive constraints on 77 
and ^ respectively, but the latter has too slow a rate to be important. 



1. Velocity dispersion of photons 

Gamma-ray bursts (GRB's) are explosive extragalactic events that release a large number of high energy photons 
with a flux that varies rapidly in time. It was therefore realized p3 . l70| that they can provide interesting constraints or 
possible observations of Planck scale suppressed Lorentz violation in the dispersion relation for photons (a possibility 
noted long ago in ) . The reason is that while propagating over such a long distance even tiny differences in group 
velocity could produce detectable time differences between the arrival at Earth of photons of different energy. 

For photons with Lorentz breaking dispersion relations of order n, ^ is related to the fractional variation in group 
velocity by 

2 M"-^ Ac 

n-1 c ■ ^^^^ 

An upper limit on the difference in arrival times of photons from the same event provides an upper limit on the relative 
speed difference, if one assumes there is no conspiracy of different emission times cancelling different propagation times. 
Together with the energies of the different photons, such observations provide a constraint on |^|. 

The strongest constraint available today comes from ORB 930131 a gamma ray burst at a distance of 260 Mpc 
that emitted gamma rays from 50 keV to 80 MeV on a timescale of milliseconds [tJ . Schaefcr 75] finds the upper 
limit Ac/c < 9.6 • 10~^^ for photons of energy ki = 78.6 MeV, and = 30 keV. This yields the constraint |^| < 122 
for n = 3. This is weaker than the constraint we have from photon annihilation, hence time of flight data do not at 
present strengthen our constraints for n = 3. For rt = 4 dispersion the bound on |f | is on the order of |^| < 10^^, so 
we get no interesting constraint for n > 3. The situation for n = 3 will be significantly improved in the future thanks 
to GLAST, the gamma ray large area space telescope, which should be able to set limits of order unity on ^ 0|. 



2. Photon Splitting 

The photon splitting processes 7 ^ 27 and 7 — > 37, etc. do not occur in standard QED. Although there are 
corresponding Feynman diagrams (the triangle and box diagrams), their amplitudes vanish. In the presence of 
Lorentz violation these processes are generally allowed when ^ > 0. However, the effectiveness of this reaction in 
providing constraints depends heavily on the decay rate. We now give an estimate of this rate, independent of the 
particular form of the Lorentz violating theory, which indicates that the rate involves at least four Lorentz violating 
factors, so is apparently too small to be relevant at observed photon energies. 

We carry out the analysis allowing for any terms in the amplitude consistent with gauge and translation invariance. 
The particular form of Lorentz violation considered in this paper also preserves rotation invariance in a preferred 
frame, however the following argument will not use that condition. Since gauge invariance is preserved, the amplitude 
for the process 7 — > N'y should arise from a term that is a scalar formed from N factors of the electromagnetic field 

(s) 

strength Fab corresponding to the external photon legs. For each photon, F^f^' ^ k^a^i,]^ where ka is the 4-momentum 
and Cb is the polarization vector. 

In the Lorentz invariant case the equations of motion imply that ka is a null vector and fcac" = 0. Energy-momentum 
conservation then implies that these 4-momenta are all parallel, so being null they are orthogonal to each other and to 
all the polarization vectors. The rate thus vanishes for two different reasons. First, since the momenta are necessarily 
all parallel, the phase space has vanishing volume. Second, the rate must be a scalar formed by contracting these 
four field strengths using only the metric. Any such contraction vanishes since it must involve contractions of the 
momenta with each other or with the polarizations. Hence the amplitude vanishes. In the case of an odd number of 



Sarkar |72| has criticize d t he use of this pa rtic ular gama ray burst since this object has no measured redshift, and hence an uncertain 
distance. Other bursts l70l or blazar flares I73l give somewhat weaker constraints. 
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photons, another reason for vansihing of the ampHtude is Furry's theorem, which states that the sum over loops with 
an odd number of electron propagators vanishes. 

If there is Lorentz violation then none of the above reasons for a vanishing rate apply. First of all the A^-odd 
amplitudes are no more guaranteed to vanish. Indeed for sufRciently general implementations of Lorentz violation 
the Furry theorem can be violated (see e.g. the discussion of the Furry theorem and its violation in the extended 
QED 77]). Secondly, the contractions of the field strengths might involve not just the metric but also a Lorentz 
violating tensor (for example u°-u^ in the rotation invariant case, where u°- is the unit timelike vector specifying 
the preferred frame.) Finally, in the presence of Lorentz violation the photon four- momenta are in general not null 
vectors hence they need not be parallel and they need not vanish upon contraction. (To satisfy energy-momentum 
conservation ^ must be positive.) 

In order for the phase space to not have vanishing volume, at least one of the 4-momenta must involve a Lorentz- 
violating factor 5 — ^(fc/Af)"^^. This is not enough for the amplitude to not vanish however. For 7 — > with 
TV = 3 or 4 the contraction of the 3 or 4 field strength tensors F^^ ~ fc[a£f)] using only the metric involves at least two 
vanishing contractions, and for larger N there are more. One of those vanishing contractions can be rendered nonzero 
by the single Lorentz violating factor already invoked on an external photon momentum, but the other one requires 
either another such factor, or a Lorentz violating tensor in the operator whose matrix element is being computed. 
Such a tensor comes with some coefficient with dimensions determined by the dimension of the operator. We also use 
the symbol 5 to indicate this sort of Lorentz-violating factor. 

The possible contributions to the amplitude will therefore be suppressed by at least two factors of 5. The rate 
goes like the square of the amplitude, hence we infer that at energies well above the electron mass the decay rate 
must behave as E6'^ or slower, where E is the initial photon energy. (There is an additional factor of if we 
consider standard QED diagrams for which each external photon leg comes with a factor of the electric charge in the 
amplitude.) 

The lifetime is therefore at least of order 5~^E~^ , which for a photon of energy 50 TeV is 10~^^5~'* seconds. Such 50 
TeV photons arrive from the Crab nebula, about 10"'^'^ seconds away, so the best constraint (i.e. if there is is no further 
small parameter such as or I/IGtt^ in the decay rate) we could possibly get on 5 from photon splitting is 5 < lO^^*'. 
For n — 2 this is not competitive with the other constraints already obtained. For higher n, each contribution arising 
from an operator of dimension greater than four will be suppressed by at least one inverse power of the scale M. For 
example, the contributions from 71 = 3 deformations to the dispersion relation will yield S ~ S^E/M. In this case 
the strongest conceivable constraint on ^ would be of order ^ < 10*, and even this is not competitive with the other 
constraints we have found. 



F. Combined Constraints 

Having completed our discussion of photon-electron processes we now turn to the determination of the global 
constraints that can be derived from the combination of all the above results. The photon splitting and the time of 
flight constraints are not as strict as those determined by the other considered interactions, at least for quadratic and 
cubic deformations, although in the future time of flight constraints may become competitive. 

1. n=2 

In the case of quadratic deviations only the difference f — 77 is constrained. The vacuum Cerenkov effect yields 
^ — f] > —10^^^, while photon decay provides the constraint (, — r/ < 10^^®. Together these confine ^ — 77 to a small 
neighborhood of zero. The photon annihilation "likelihood region" would just impose £, ~ rj < 10~^^, which does not 
further strengthen the constraint. 

2. n=3 

Putting together the constraints from the three photon-electron interactions previously considered we obtain a 
remarkably small allowed region in the 77-^ plane (see Figure ISJ. The photon decay and Cerenkov constraints exclude 
the horizontally and vertically filled regions respectively. The allowed region lies in the lower left quadrant, except 
for an exceedingly small sliver near the origin with < 77 ^ 10~^ and a small triangular region (—0.16 ^ 77 < 0, 
< ^ < 0.08) in the upper left quadrant. The discussion of the photon annihilation threshold in subsection IIII U j| 
indicates that, although no firm constraint can be given at present, the allowed region cannot lie too far from the 
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FIG. 8: Combined constraints on the photon and electron parameters, for the case n = 3. The regions excluded by the photon 
decay and Cerenkov constraints are lined horizontally in blue and vertically in red respectively. The region between the two 
diagonal green lines is where the threshold for the annihilation of a gamma ray with a 25 meV photon ranges from its standard 
value (upper diagonal green line) to not more than twice that value. The shaded patch is the part of the allowed region that 
falls between these gamma annihilation thresholds. The dashed line is ^ = 77. 



corridor between the two roughly parallel diagonal lines. These lines indicate where the threshold for the annihilation 
of a gamma ray with a 25 meV photon ranges from its standard value (upper diagonal green line) to not more than 
twice that value. 

If future observations of the blazar fluxes and the IR background yield agreement with standard Lorentz invariant 
kinematics, the region allowed by the photon annihilation constraint will be squeezed toward the upper line (fcth ~ ^s)- 

Time of flight constraints for high energy photons currently constrain ^ to be less than ^ 100 at best, but future 
observations should allow such constraints to further narrow the allowed region towards the origin. 

3. n=4 

The case of quartic deviations is unfortunately just mildly constrained from the available observations. The order of 
magnitude allowed for the parameters is as small as 10^^ (from Cerenkov) for the electron-photon vertex interactions. 

IV. INTERACTIONS WITH PROTONS, NEUTRINOS, AND MUONS 

We have focused so far on effects involving just electrons and photons, in order to determine the strongest available 
combined constraints. We now briefly discuss some other interactions that are realizable with a violation of Lorentz 
invariance, and which can now or in the future provide further constraints or observations of Lorentz violation. 

A. Alternative vacuum Cerenkov effects: protons, neutrinos and muons 

The former discussion of the vacuum Cerenkov effect can be applied also for any other particle that couples to 
photons, using the same kinematic equations. Since the strength of the observational constraint is determined by the 
smallness of the ratio rii? /p^^ax, smaller masses or larger energies generally lead to stronger constraints. However, in 
the case of neutral particles that couple to photons only through higher multipole moments the rate must also be 
considered. We summarize in Table Ull the values of the quantity vr? jVmax- 
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TABLE II: Typical values for different particles for the actual or potential constraints from absence of the vacuum Cerenkov 
effect. 





V 


e 




p+ 


m 


^ 1 cV 


0.511 MeV 


105 MeV 


938 MeV 


Pmax 


~ 1 TeV - 10^° eV " 


~ 100 TeV ' 


~ 1 PeV 


~ 5 ■ 10'^ eV 




n = 2 


~ 10-2* - 10-*° 


n = 2 


~3 - 10-^'' 


n = 2 


~ 10-" 


n = 2 


~ 4 ■ 10-22 


n = 3 


~ 10-* - 10-^2 


n = 3 


~ 3 ■ 10-^ 


n = 3 


~ 10-^ 


n = 3 


~ 8 ■ 10-" 


n = 4 


~10« -10-2* 


n = 4 


~ 3 • 10" 


n = 4 


~ 10^2 


n = 4 


~ 2 • 10-^ 



"Lower value is AMANDA data; largest value is potentially observable UHE neutrinos. 

'Energy expected for electrons responsible for the creation of ~ 50 TeV gamma rays via inverse Compton scattering |2^|^. 
■^Expected energies to be detected for muons produced by cosmic neutrinos. 
"^Detected in UHECR. 



1. Protons 

Very strong constraints can be obtained using the ultra high energy protons in cosmic rays, up to the GZK cutoff 
of 5 • lO'^^ eV. The identity of these particles has been called into question by the candidate events beyond the GZK 
cutoff as described in Sect. IIVBI However, even if the highest energy events do not originate with protons, there is 
strong evidence that protons up to the GZK cutoff do exist in cosmic rays p^.-'-^ 

The rate of vacuum Cerenkov radiation from charged particles is irrelevant for the determination of constraints 
since it is very high. (See Sect. IIII B 21 ) For the parameter region where the threshold occurs with emission of a zero 
energy photon, the proton can presumably be treated as a point charge so the threshold relations previously obtained 
for electrons are directly applicable using the proton mass in place of the electron mass, and the parameter rjp from 
the proton dispersion relation in place of rjf,. This region of parameter space is described in section Fill Bl 

For parameters where a hard photon is emitted at threshold, the role of the partonic structure of the proton needs 
to be examined, which we have not done. It may turn out that the threshold can be determined by the quark 
dispersion relation rather than that of the proton. If so, it would be the quark deformation parameter rjq rather than 
•qp that is constrained by observations of non-decaying high energy protons, and one would need to use the quark mass 
and energy in the threshold relations. In this case the proton may be destroyed rather than just slowed by vacuum 
Cerenkov radiation, however that distinction is irrelevant for the determination of constraints, since either way high 
energy protons would not travel long distances. 

In estimating constraints we ignore here the possible role of partonic structure, and simply use the proton mass 
and energy in the threshold formulae derived in section FlII BI for point particles, with the understanding that for hard 
emission thresholds the constrained parameter may be rjq rather than rjp, and the numbers may be off by a few orders 
of magnitude since the quark mass and energy were not used. 

Using the GZK cutoff (5 • 10^^ eV) for the highest energy protons we obtain the following constraints relating the 
parameter ^ in the photon dispersion relation and rjp in the proton dispersion relation. For a quadratic deformation 
of the dispersion relation (n = 2) the bound is ry^ — ^ < 4 • 10"^^. For cubic deformations (n = 3) the constraints on 
parameter space have the same form as represented in Figure In the case of the proton the quantity rn^/p^^^, is 
of order lO-^"* compared with IQ-'^ in the case of 100 TeV electrons, which means that the boundaries of the allowed 
region are closer to the ^ axis in the upper half plane and to the diagonal in the lower half plane. However, the 
qualitative nature of the allowed region is identical. A good constraint is even obtained for the case of quartic (n = 4) 
deviations. As shown in section IlII B II it is the quantity rn^/ [p^^^ that determines the strength of the constraint 
in this case. For 5 • 10^^ eV protons this is approximately 10"^, still much less than unity and a much better figure 
than the 10^^ obtained for the 100 TeV electron. For n — b deviations the strength of the constraint is determined 
by 'm'i / {p%^^y-) ~ 10'^, hence one does not obtain even order unity constraints on the coefficients. 



Note added in proof. A recent analysis 85j argues that there are insufficient statistics to establish the GZK cutoff at this time, hence 
the existence of these protons cannot yet be regarded as established. 
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2. Neutrinos 

In the standard model the vacuum Cerenkov reaction with neutrinos, — > + 7, is not allowed due to energy- 
momentum conservation - whether or not the neutrinos are massive. If they are massive energy-momentum conser- 
vation cannot be satisfied at all. If they are massless it can only be satisfied if all three particles are strictly parallel, 
yielding no phase space for the reaction. (Since there is good evidence that neutrinos have mass, we will assume this 
for the rest of the discussion.) Energy- momentum conservation is the only obstruction for this reaction, since although 
the neutrino is neutral there is a nonzero matrix element for the process. In particular there are two channels: the 
charge radius interaction and, if massive, a magnetic moment interaction (see e.g. HJ])- We therefore see that, as for 
charged leptons, Lorentz violating dispersion relations can allow the reaction to happen. 

In order for the neutrino Cerenkov reaction to give strong constraints on Lorentz violation two conditions must be 
satisfied: (1) the energies where Lorentz violating terms are comparable to the neutrino mass term in the dispersion 
relation must be accessible to observation, and (2) the rate of the reaction must be high enough so that it would 
significantly affect the propagation of observed neutrinos. The first condition is already met since the relevant energy 
where Lorentz violation becomes important is 100 MeV (see TableHl, while Super Kamiokande has detected neutrinos 
over 100 GeV |73] and the AMANDA detector has seen neutrinos up to a few TeV j79j. The second condition is more 
problematic since both the charge radius and magnetic moment channels are very strongly suppressed. 

The best case for current observations would be using AMANDA, since the neutrinos have the highest energy 
and travel the diameter of the earth after being produced in the atmosphere above the North Pole. We have not 
carried out a detailed analysis, but an estimate given below suggests that the Cerenkov rate is not high enough to 
produce interesting constraints with these neutrinos. The energy loss rate depends strongly on the energy however, 
so atmospheric PeV neutrinos, which are likely to be detected by AMANDA or IceCube (seee.g. jSm), may provide 
constraints. (The same experiments should detect PeV muons as secondary products of the neutrinos, which would 
also provide an interesting constraint as seen in Table UTI) Still higher energy neutrinos, up to perhaps 10^*^ eV, are 
expected either as cosmic ray primaries or as a byproduct of cosmic rays 10, 81]. Such high energy neutrinos could be 
detected by AMANDA 82] , and they could be observed via horizontal or possibly upward air showers using existing 
detectors like HiRes or future ones such as the Telescope Array [s^, 0] . 

In Table we summarize the typical constraints one can expect from neutrinos in the above mentioned range of 
energies provided the rate is high enough. Remarkably, the combination of high energies and low mass could give for 
cosmological neutrinos [E^ 3> 1 PeV) stringent constraints {r^u <C 1) for deviations up to n = 6. 

A calculation of the neutrino vacuum Cerenkov rate is beyond the scope of this article but we provide here a rough 
estimate that may provide some guidance in this problem. We saw before (Section IIII Bp there are two types of 
Cerenkov thresholds depending on the values of ^ and 77^: the "soft" one which occurs with emission of a zero energy 
photon, and the "hard" one in which a photon with energy comparable to the incoming particle is emitted. The decay 
rate will be much greater in the hard threshold case, so we consider that here. (The soft threshold case may still be 
relevant well above threshold.) 

To show that the rate might be fast enough to provide a useful constraint it suffices to examine the charge radius 
interaction. This occurs via the emission of a virtual W-boson, hence the amplitude goes like Kfc|/M^, where k is a 
small numerical factor (~ 10^®) coming from coupling constants and integration measure, k\ is the square of the photon 
four-momentum, and Mw is the W-boson mass. We thus estimate the rate for Cerenkov emission from a neutrino of 
very high energy £^ to be F ~ [nkH M'^)'^ E^. (The factor of E^ is determined by the phase space integration, which 
does not involve any Lorentz violating factors well above threshold.) With Lorentz violating dispersion of order n we 
have k\ = ^/c", hence the rate goes like F ~ {K£,k" /M^)'^E^. Taking the photon energy to be of the same order as 
the neutrino energy k ~ E,^, this gives a lifetime for emission r ~ ^^^(i?iy/PeV)~'^"+^-' x io26"-86 seconds. If correct 
this would be short enough to yield interesting constraints for n — 3 using atmospheric PeV neutrinos travelling 
through the earth, since their transit time is of order 10~^ s. 

As a final remark, we note that the related process of photon decay to two neutrinos could also take place in the 
presence of Lorentz violation. This would yield strong constraints on f and rj^ provided the rate is high enough. The 
above estimate suggests that for multi-TeV photons from cosmological sources the rate would indeed be high enough 
for n = 2, 3. 



Note that we cannot constrain ^ as much as the threshold would indicate, since for extremely small ^ the decay rate eventually gets too 
small. 
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B. The GZK cutoff 



The presence of the GZK cutoff on the ultra high energy (UHE) proton spectrum is due to pion photoproduction: 
7P — > p7r°, as previously discussed in the Introduction. The observation of this cutoff also gives constraints on the 
Lorentz violating coefhcients. Current data from the HiRes, Fly's Eye and Yakutsk experiments strongly indicate that 
the GZK cutoff is present at a cosmic ray energy of 5 • 10^^ eV While AGASA reports a number of extra events 
beyond the expected flux of high energy cosmic rays above 10^" eV, below 10^" eV AGASA also shows evidence for 
the GZK cutoff (see e.g. Figure 1 of [23|)- Unfortunately, the experimental data are strongly affected by the uncertain 
energy calibration of each experiment. A systematic analysis of the data allowing for various calibrations is outside 
the scope of this work, so for now we assume that the published energy calibrations are correct. 

We constrain rjp^rjT^ by determining where the induced modification of the cutoff would disagree with the data. 
(The incoming photon has low energy and so no useful constraints on ^ are obtained.) The 77^, rj-n constraints are 
quite strong (on the order of 10~^° for n = 3) due to the high energy of the reaction. 

In the standard Lorentz invariant theory the threshold energy for pion production is £'th = mT^{2mp + r7i.n.)/4a;o, so 
a photon with energy ~ 1-3 meV is at threshold with the proton at the GZK energy. In order to give a constraint on 
Lorentz violations we consider raising or lowering the UHE proton at threshold with the same luq. This is equivalent 
to changing the GZK cut-off as we are modifying the UHE proton energy that interacts with the relevant CMBR 
photons responsible for the Lorentz invariant GZK effect. 

Examination of the data plot in Fig. 2 of Ref. [23| reveals that if the cutoff were shifted via a Lorentz violating effect 
down to 2 • 10^^ eV or up to 7 • 10^^ eV then the theoretical predictions would no longer agree with the data at above 
a 2a confidence level. This energy range therefore provides constraints on rjp^rj-^. From the threshold theorems of 
we again know that in the threshold configuration where the GZK reaction begins to occur the incoming proton and 
photon collide head on and the outgoing proton and pion 3-momenta are parallel. Energy-momentum conservation 
in this configuration and the dispersion relations give an equation similar to equation H46|) for photon annihilation. 
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where x = q/p, and p and q are the initial and final proton 3-momenta. 

For n = 2 the threshold analysis has already been done by Coleman and Glashow 17] leading to a constraint 
tItt — r/p < 5 X 10^^'* [^/'^]^ for a target photon uj, where uj = fcTcMB = 0.235 meV. 

For n = 3, the presence of the pion in this equation complicates the analysis as there is an additional mass term 
and the fi nal p articles are not interchangeable. The case of equal coefficients (rjp = 77^) has been studied analytically 
in |18L I19I I20I l2a. |27| and numerically in 29] . Here we numerically find the thresholds for the GZK reaction allowing 
for unequal coefficients in the case n — 3. The region in the r]p, rj-^ plane where the thresholds arc in the allowed range 
discussed above are shown in Figure El in which the axes are in multiples of 10"^". 
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FIG. 9: The range of rip,rjn for n — 3 dispersion modifications where the GZK cutofi^ is between 2 ■ 10^^ eV and 7 ■ 10^^ eV. rjp 
and r]^ are in multiples of 10~^°. 
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We turn now to the question of the extra AGASA events above 10^° eV. The AGASA data is sparse in this 
energy range, and there is not a large, precise data set from other experiments with which AGASA disagrees. The 
uncertainties in all experiments are large enough that a modified theoretical theoretical spectrum could possibly agree 
with all experiments at the la level. One cannot therefore simply disregard the possibility that the flux above 10^° 
eV is in fact higher than the standard theoretical prediction. 

Previous authors have suggested that the AGASA events above the GZK cutoff could be explained by an upward 
shift of the GZK cutoff induced by Lorentz violation El El El IM El El , however this seems incompatible with 
current data since the cutoff is seen. Another, more subtle possibility is that these events are related to the existence 
of an upper threshold. We have checked numerically that no upper threshold exists below 10^° eV within the allowed 
region of Fig.|51 Nevertheless, the phase space for a reaction begins to close up before the upper threshold is reached. 
The reduction in phase space would in turn reduce the rate of the GZK reaction leading to a higher than expected count 
of events at high energies. If the lower threshold were dramatically modified whenever there is an upper threshold 
then this scenario could not explain the data. However, this is not the case - there are choices of r/p, t^tt such that an 
upper threshold exists and the lower threshold is only slightly modified. Since the lower threshold modifications can 
be small, the experimental signature of the GZK cutoff could remain unchanged near 5 • 10^^ eV while the intensity 
of the spectrum at high energies increased from its Lorentz invariant prediction. This scenario could perhaps explain 
the AGASA data and be compatible with other experiments if the upper threshold is low enough that there is a 
significant phase space reduction just above lO^*' eV. The range of Typ,??^ for which this effect could occur is given in 
Figure El There we consider the range of parameters for which the lower threshold still lies between 2 • 10^^ eV and 
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FIG. 10: The range of r/p, r;^ for n = 3 dispersion modifications where the GZK cutoff is between 2 • 10^^ eV and 7- 10^^ eV and 
the upper threshold exists below 10'^^ eV. ??p,7?7r are in multiples of 10~^°. Note that this region does not include the origin, as 
there is no upper threshold in the Lorentz invariant case. 

7 • 10^^ eV but an upper threshold exists below 10^^ eV (so that the induced reduction in phase space could affect the 
AGASA data). 

We have not considered the constraints that can be obtained from the GZK cutoff in the case n — 4 but it is clear 
that they are interesting in this case as well, since 1/pmax = Af/Pmax = 2 • 10®. The n = 3 constraints are of order 
10"^*^, hence the n = 4 ones are of order 10^^. 

C. Neutron stability— Proton instability 

If there are different dispersion relations for protons, neutrons, positrons, and neutrinos then protons may be 
unstable and decay to neutrons at sufhciently high energies. For n = 2, Coleman and Glashow E3 have shown this 
explicitly. If n > 2 then the analysis becomes more complicated, but there exist parameters for which the neutron is 
stable and the proton decays. For example, consider a neutron with rjn — ^1 and n — 3 dispersion relation, and an 
unmodified proton, electron, and neutrino. At momenta p„ > mf/^ the neutron energy-momentum vector becomes 
spacelike. Since the energy momentum vectors of the other particles are still timelike it is impossible to satisfy energy 
momentum conservation hence neutron decay does not occur above this energy. This opens the possibility that ultra 




27 



high energy cosmic rays are neutrons rather than protons, in which case the GZK cutoff for these cosmic rays is 
irrelevant since neutrons do not interact strongly with CMBR photons. The presence of the observed cutoff can thus 
be used to constrain the parameters further, and one could also contemplate the possibility that the AGASA events 
above the GZK cutoff are present because the neutron becomes stable just above the cutoff. 

V. COMBINED CONSTRAINTS FOR A UNIVERSAL LORENTZ BREAKING DISPERSION 

RELATION 

Although we have considered so far the case of different Lorentz violating parameters rja for different particles, it 
may be that the underlying quantum gravity physics selects a universal deformation parameter rj. This indeed was 
the ansatz considered in most of the previous literature. We therefore consider now this special case with n — 3 
deformations. 

We start by considering the photon-electron interactions. From Figure ^ ^^nd equation (|23|l we see that the 
Cerenkov effect limits the available values on the diagonal to a semi-infinite line with 77 < m? /2p^^^. This corresponds 
numerically to 77 < lO"^** if we consider the constraint provided by the observation of ultra-high energy protons in 
cosmic rays. The analysis of photon decay shows that the permitted values on the diagonal ^ = rj are restricted to 
the semi- infinite line 77 < Sm^/fej^^^^ ^ 8 • 10~^ (using the observation of 50 TeV gamma rays from the Crab nebula). 
The observation of photon annihilation provides, as previously discussed, a more uncertain constraint. Nevertheless 
for our purposes it is enough to take into account that some absorption is detected for gamma rays at least up to 10 
TeV. We can then take as a definite constraint the line for the existence of a lower threshold as shown in Figure |S1 
This line meets the diagonal at 77 = —32/27. The problem is now to decide for which ujq we are sufficiently confident 
the photon annihilation still takes place. As a reference value we take here again the = 25 meV photon previously 
considered. In this case the region of existence of a lower threshold for the photon annihilation limits the value of i] 
to the semi-infinite range 77 > -2.3 x (32/27) « -2.7. 

If the GZK cutoff is confirmed that would establish with certainty that at least some of the UHE cosmic rays are 
indeed protons. Moreover it would also provide a correspondingly strong constraint on negative values of 77. If the 
GZK cutoff is within order unity of it's Lorentz invariant value, 77 is constrained to be |7;| < 10^^^. (Note that this 
constraint on 77 is so strong as to exclude the region of upper threshold for the GZK process shown in Figure [TUl 'l 
The upper bound might be further pushed toward zero if one takes into account the Cerenkov effect of high energy 
neutrinos (see Table Hlj) . 

VI. DISCUSSION 

In this paper we have performed a systematic analysis of the effects of Lorentz violating dispersion on particle 
reactions, allowing for unequal deformation parameters for different particles. We have analyzed the threshold kine- 
matics and combined the observational constraints where possible. Even when suppressed by the inverse Planck mass, 
such Lorentz violation can lead to radically new behavior in the kinematics of particle interactions at much lower 
energies. Reactions previously forbidden can be allowed, lower thresholds can be shifted and upper thresholds can be 
introduced. The presence of upper thresholds is a feature of Lorentz breaking physics that is not present in Lorentz 
invariant physics and which can be relevant for observational constraints.^'* Furthermore, we have found that for 
interactions with identical final particles, the final momenta can be distributed asymmetrically at threshold. While 
this is a straightforward consequence of the kinematics, it has been previously overlooked in the literature, probably 
because it is alien to Lorentz invariant physics. 

Using these kinematical results, we have seen that a conservative interpretation of observations puts strong con- 
straints on the coefficients 77 and ^ of order E/Mp modifications to the electron and photon dispersion relations. The 
allowed region includes ^ = 77 = — 1, which has been a focus of previous work pol l25l l2^. The negative quadrant 
has most of the allowed parameter range. Note that in this quadrant all group velocities are less than the low energy 
speed of light. For modifications of order (E/Mp)'^ there are no significant constraints in the electron-photon sector 
derivable from current observations, due to the fact that the energies of observed particles are too low. However 
reactions such as proton Cerenkov (in vacuo) or pion production by cosmic rays, for which we have data at much 
higher energies, can provide good constraints for (E/Mp)'^ modifications (although for different particle deformation 
parameters). Ultra high energy cosmological neutrinos may also provide good Cerenkov constraints at this or even 



More complex dispersion relations can lead to multiple thresholds 16^ with could have further observational effects. 
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higher orders, since the neutrino mass is much smaller than that of any other particle. The interaction amplitudes 
are very suppressed however, so it is necessary to accurately calculate the rate and compare it with the travel time of 
the neutrino. 

There are a number of ways to improve the constraints on 0{E/Mp) modifications from electron-photon interactions. 
Higher energy electrons would not help much since the Cerenkov constraint is already strong, while finding higher 
energy undecayed photons would squeeze the allowed region onto the line ^ = 77 of Figure |H1 To further shrink the 
allowed segment of this line would require improved knowledge of the infrared background and a reconstruction of the 
source spectrum from the observed gamma rays in the presence of Lorentz violation. Also, the constraint from time 
of fiight measurement may become competitive using improved detectors. 

Other constraints may be provided by additional interactions not considered here. For example, a possible upper 
threshold for e+ ^ 27 cannot provide a competitive constraint in astrophysical observations since there are other 
processes by which observed high energy photons can be produced. However, if future electron accelerators can reach 
energies above 10 Tev then one can expect to get a good constraint from this reaction. In addition there may be other 
reactions for which upper thresholds can produce useful constraints at or near currently observed energies. Reactions 
involving more than two types of particles, such as z/ — > e~ , could also give constraints. It may be possible that 
by considering a number of such reactions a multi-dimensional parameter space can be usefully constrained. 

The idea motivating our work is that Lorentz violation may be a consequence of quantum gravity, in which case 
the natural scale for the Lorentz violation is the Planck scale. If, as in braneworld scenarios, the quantum gravity 
scale were to be around a TeV, then the natural scale for Lorentz violation induced by quantum gravity would be the 
TeV scale. Clearly, the only way such Lorentz violation could be compatible with observations is if it were extremely 
suppressed compared with this natural scale. This suggests that either TeV scale quantum gravity is wrong, or it 
does not violate 4d-Lorentz invariance. 

In conclusion, the absence of anomalous observations provides stringent constraints on the possibility of Lorentz 
violation originating at the Planck scale. This in turn gives important information as to the viability of quantum 
gravity theories that predict 4-d Lorentz violation. We can expect that, as better data at higher energies becomes 
available, even stronger constraints will be imposed or, alternatively, positive signatures of Lorentz violation may be 
found. Either way, it is clear that a useful tool for the phenomenological investigation of quantum gravity is now at 
hand. 



Acknowledgments 

We are grateful to G. Amclino-Camelia, A. Celotti, R. Mohapatra, I.Z. Rothstein, and F. Stecker for useful dis- 
cussions. This research was supported in part by the National Science Foundation under Grant No. 9800967 at the 
University of Maryland, Grant No. PHY99-07949 at the KITP, and the Erwin Schrodinger Institute. 



APPENDIX A: PHOTON ANNIHILATION THRESHOLDS 



In this appendix we work out the lower and upper thresholds for the process 7 7 ^ e as a function of the 
Lorentz- violating parameters r/ and ^. Our starting point is the kinematic equation (|4(j|l derived in section Fill Dl 

= F{k, y) + (^e + ^1^) 2/ (1 - jy) - w (1 - y) [i''-'^ + (1 - yt-'''] • (Ai) 

Here m is the electron mass, k is the magnitude of the incoming hard photon momentum, loq is the soft photon 
energy, y = p/k, where p is the magnitude of the electron (or positron) momentum, and the threshold configuration 
of antiparallel incoming photons and parallel outgoing electron-positron pair has been imposed. This equation follows 
from (i) energy-momentum conservation, (ii) the dispersion relations for the particles, and (iii) the threshold configu- 
ration. To find the lower and upper threshold for given values of 77 and ^ we must determine the minimal or maximal 
k for which the reaction can occur. According to the threshold theorem (c/. Sect. IIII A|) these k values always occur 
with what we just called the threshold configuration, hence we must determine the minimal or maximal k for which 
there is a solution (fc, j/) to the kinematic equation IjAip with y in the range [0, 1]. In the Lorentz invariant case the 
threshold always occurs with the symmetric configuration y = 1/2, however in the Lorentz- violating case this is not 
always true. 

In order to derive results applicable to any value of the soft photon energy ljq and "electron" mass m, we introduce 
scaled variables 
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where /cli is the standard lower threshold m? /ujq. In terms of these scaled variables the equation m '^k"F{k, y) = 
takes the form 

G{P, y) ^ a.n{y) (3" +j{y) (3 -1 = 0, (A3) 

where 

a^iy) = y{l - y) (j - v [y^-' + (1 - (A4) 

and 

7(y) = 4 2/(l-y). (A5) 

Equation (jA3p is a generalization of those derived by Aloisio et al. [2^ for the specific cases of n = 3,4 with purely 
symmetric configurations {y = 1/2) and equal deformation coefficients = 

Figure ITTI shows the general behavior of (?(/?, y) (jA3|l for any fixed y (and therefore fixed 7(1/)) and different values 
of a. From the plot we see that for any y there are either one, two, or zero solutions to the kinematic equation for 




FIG. 11: Behavior of G{f3,y) JA3II for fixed y and different values of alpha. 



j3. If there are two solutions, of course only the lower one is a candidate for the lower threshold. The upper one can 
correspond to an upper threshold, that is, to the highest available value of (3 for which the reaction is kinematically 
allowed. To our knowledge, the possibility of upper thresholds has been overlooked in all of the previous literature 
except for Kluzniak (^^e discussion at end of section IlII D 3p . 

For each y there is a maximal value of the lower root /3 which occurs when a„ is such that the curve described by 
Eq. l|A3p is tangent to the j3 axis. This occurs for 

a'r^ = -r- , (A6) 

and, at this tangency point, 

/3tang=7'';^- (A7) 

If /3 is to be a lower threshold it must lie below this tangency point. 

Given values for ry and ^, Eq. I A3I implicitlv defines zero, one, or two (real positive) solutions for /3 as a function of 
y. A lower threshold corresponds to the global minimum of (3. Our strategy is to first find the local minima and, if 
there is more than one, determine which is the global minimum. 

A local minimum of (3{y) is characterized by = and (3yy > 0, where the subscript y denotes derivative with 
respect to y. To find the corresponding conditions on G we use the fact that G{(3{y),y) is identically zero, hence its 
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total derivative with respect to y vanishes: G/^Py + Gy = 0, where the subscripts on G denote partial derivatives. The 
second total derivative of G{(3{y), y) with respect to y also vanishes. At a stationary point where (3y = 0, this yields 
Gf3l3yy + Gyy = 0. Thus thc conditions for /?(?/) to be local minimum are 

f3y^-Gy/Gp^0 and Pyy ^ ~Gyy/Gp > 0. (A8) 

It is clear from Fig. ^]that is always positive at the lower root of G = 0, and vanishes only at the tangency value 
of a. Since only the lower root can be a lower threshold, we have two necessary conditions to be a lower threshold: 

Gfj>0 and Gyy <0. (A9) 

We now use these considerations to find the thresholds for n — 2,3. 



1. Photon annihilation thresholds for n — 2 

For n — 2 Eq. (|A3|) reduces to 

G = y{\- y) [{I- v)f3' + 4/?] - 1 = (AlO) 
There is only one extremum, at y = 1/2. Substituting y = 1/2 into Eq. (jA10|) yields 

^=?, + A^. (All) 

The ^-intercept, = 4(1 — /3)//3^, decreases monotonically for f3 < Aang and increases monotonically for (3 > Ptang- 
d^o/df3={f3-2)(4/f3^). 

The values of /3 less than /3tang = 2 are candidates for a lower threshold. The contours I|A11|I for /3 < 2 are parallel 
straight lines of unit slope, whose ^-intercept goes monotonically from oo to — 1 as /3 goes from to 2. Since these 
lines do not cross there is only one candidate threshold for each pair [rj, ^) , hence these lines indeed give the contours 
of the lower threshold. Thc reaction is forbidden below the line ^ = t^ — 1. Thc highest lower threshold is given by 
k = 2ku — 2m^/cjo- 

The values of P greater than /3tang — 2 are candidates for upper thresholds since they arc local maxima of (3 among 
the threshold configurations. The contours for these are also given by the straight lines ljAll|l . with ^-intercept that 
goes monotonically from —1 to as /3 goes from 2 to oo. No other candidate for the upper threshold exists at a 
given value of f , rj, so to check whether these contours actually represent upper thresholds we need only verify that 
an upper threshold exists at all. That is, we must rule out the possibility that there are configurations for which the 
annihilation process occurs with arbitrarily large incoming photon energy. To do this we examine the limit of large 
k. Then the soft photon is irrelevant, and the question is the same as whether an arbitrarily high energy photon can 
decay to an electron-positron pair. We already determined in the section on photon decay that it occurs only above 
the diagonal, i.e. for ^ > rj. Hence the process does not occur at arbirarily large energy below the diagonal, which is 
where all the candidate upper thresholds lie. Thus these candidates are indeed all upper thresholds. 



2. Photon annihilation thresholds for n — 3 



Since the two final particles are interchangeable (|A3(I is symmetric about y — 1/2. For higher n we can reduce the 
order of this equation by introducing the variable z = (2y — 1)^ that also has this symmetry, as we did for the case 
of photon decay in section IlII CI The physically relevant range of z is < 2; < 1. With this change of variables for 
n — 3 (|A3f) can be written as 



G 



0- 



'3 r 



(e + 4//?2)(l-z)-(?7/2)(l-z2) 



1 = 0. 



(A12) 



There are now two extrema, one at z = which corresponds to the usual symmetric y = 1/2 case, and the other 
where G, = 0, i.e. at 



C + 4//32 



(A13) 
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which corresponds to an asymmetric configuration in which the outgoing particles have different momenta. 
The solution of (|A12I) for ^ in the symmetric case (z — 0) with (3 = (3s yields 

The ^-intercept, — 4(1 — (3)/P^, decreases monotonically for f} < /3tang = 1.5 and increases monotonically for 
/? > Aang d^o/df3 — {(3— 1.5)(8//?^). Only values of /3 less than the tangency value 1.5 are candidates for a symmetric 
lower threshold, while values greater than 1.5 are candidates for a symmetric upper threshold. 

The symmetric case can only be a lower or upper threshold when the inequality ^ > holds. We can see this 

by imposing the conditions for a local minimum or local maximum respectively: f3yy > when /3 < 1.5 or Pyy < 
when P > 1.5. Since G/3 > in the first case and G/3 < in the second case (as can be seen from Fig.^fl, (|A9|) shows 
that both cases require Gyy < 0. To evaluate Gyy we note that 

(f/dy'^ = 16z{(f/dz^) + 8{d/dz), (A15) 

hence Ci"^ = SCi"^ = -8/3(1 + 63^4), which is negative only if f > -4/^^ 
The solution of (|A12p in the asymmetric case {z = Zg) with [3 — Pa yields 



The asymmetric case only exists when ^ < 0, and only the positive square root is physically relevant, since 
gives = 1 — \f^8/rf(3^ and z must be less than unity. Also, Za must be positive, so the asymmetric case is only 
relevant when ^ < — 4//3^. 

The asymmetric stationary point has Gy'^y — IQzaG^zz = ^Zarj0^ < 0, hence it is a local minimum if and only if 
Gp > 0. This corresponds to the inequality ^ > rj + 4/3(3'^ . It is a local maximum when the opposite inequality holds. 
For (3 < 1.5 the asymmetric curve represents a local minimum everywhere since it is above ^ = rj + 4/3/3^ everywhere 
in the physical region ^ < — 4//3^. For /? > 1.5 the asymmetric curve crosses below ^ = rj + 4/3/3^ while still in the 
physical region, where it represents a local maximum. 

The symmetric line ljA14|l and asymmetric curve (IAI6II meet at (77,^) = (— 8//3'^, — 4//3^) and are tangent there, as 
shown in Figure IT21 Above this meeting point only the symmetric solution is a candidate threshold, and below this 




FIG. 12: The (straight) symmetric (I A 1411 and (curved) asymmetric 1IAI6I 1 contours for a fixed value of (3. The unphysical part 
of the asymmetric curve is dotted, and the part of the symmetric line that is not a local maximum or minimum of l3 is dashed. 
The "/3-curve" is the two joined solid segments of these contours, and indcates points where /3 is a candidate for the lower or 
upper threshold. The joining point is (rj,^) = (— 8//3^, — 4//3^). 
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point only the asymmetric solution is. As (3 varies, the curve traced out by these meeting points is given by 

§oin - -i-v)'/'- (A17) 

We shall use the name "/3-curve" for the joined curve that is the symmetric line above and the asymmetric curve 
below ^= -4/^2 _ 

a. Lower threshold for photon annihilation, n = 3 

To find the contours of constant lower threshold in the rj~£^ plane we proceed as follows. First we choose a value of 
(3 < 1.5, and consider the corresponding /5-curve. The points on this curve are the only candidates for the threshold 
to be p. To determine if the threshold actually is /3 at a given point we must determine whether or not there is 
a solution to (jA12|l with a smaller value /3< at the same point. In other words, we must determine if a /3<-curve 
could cross the original /3-curve. In fact it cannot. The /3<-curve starts out above the /3-curve at ^ = (since 
(1 — /3<)//3< > (1 — P)/ (i^)- In the symmetric sections the slopes are both equal to 1/2. In the asymmetric section 
the slope computed from ljA16|) is d^/dp — 1 ~ 4(— Sry/J'^)"^/^, which is always greater than 1/2 in the region below 
^ — — 4//3^ and is greater for larger f3 at fixed rj. Hence the /3<-curve is everywhere above the /3-curve, so the curves 
never cross. The /3-curves with (3 < 1.5 thus give the lower threshold. 

Now assume that (3 is greater than 1.5, so that only the the asymmetric configuration can be a lower threshold. In 
this case there are /3<-curves with smaller values of /3< that cross the /3-curve. In particular, the asymmetric part of 
the /3-curve crosses the symmetric /3 = 1.5 line from below and then goes on to cross lines of yet smaller f3 above that 
before leaving the local minimum region. Thus /3 cannot be the global minimum above the symmetric f3 = 1.5 line, 
even though (3 remains a local minimum up to when it crosses below the line £. =7} + 4/313'^. 

The only remaining question is whether the different asymmetric /3-curves for (3 > 1.5 can cross below the symmetric 
(3 = 1.5 line. In fact they cannot. It can be shown that the terminus on the (3 = 1.5 line moves to larger values of rj 
as f3 goes from 1.5 to oo. (See Figure 0) Since the slope of the asymmetric /3-curves increases with /3 (as discussed 
two paragraphs above) in the region below the symmetric (3 = 1.5 line, the curves for different (3 do not cross. 

In summary, we have shown that the regions where the symmetric and asymmetric lower thresholds exist take the 
form shown in Figure 2| The contour lines of constant threshold are the straight lines (|A16p of slope 1 /2 in the 
symmetric region, and are given by ljA16|l in the asymmetric region. These contours are shown in Figure |S1 

b. Upper threshold for photon annihilation, n = 3 

We now turn to analysis of the upper thresholds. Our first step is to ascertain in which region of the 77-^ plane 
an upper threshold exists. As discussed in the n = 2 section, this can be done by examining the limit of large fc, in 
which the annihilation process becomes indistinguishable from photon decay. The decay process is forbidden below 
the broken line given by ^ = v/2 for = rj > and ^ = rj for 77 < 0. Above this broken line there is a lower threshold 
and no upper threshold. Thus an upper threshold exists below this broken line anywhere a lower threshold exists. 

The candidates for upper threshold contours are the sections of /3-curves with /3 > 1.5 that satisfy the conditions 
for being a local maximum. On the symmetric segment this imposes no restriction, but on the asymmetric segment it 
requires that the curve lie below the fine £, — rj + 4/3/3^. On the other hand, we just argued that an upper threshold 
exists only below the diagonal for negative rj, which is a more restrictive condition. Moreover, it can be checked that 
these sections of the /3-curves do not cross anywhere in the region of upper thresholds, hence they are indeed the upper 
threshold contours in that region. Only for /3 > 2 does the asymmetric section have a segment below the diagonal 
before leaving the physical region ^ < — 4//3^, as illustrated in Fig. ^1 The regions of symmetric and asymmetric 
upper thresholds thus take the form shown in Figure El The the boundary of the lens shaped region next to the 
diagonal is determined by the curve (jA17|l consisting of the points where the symmetric and asymmetric segements 



This analysis was carried out using the truncation discussed in sectior lTll Al At very high energies or at very large values of the Lorentz 
violating parameters there will be deviations from this behavior, but for practical purposes these deviations are not relevant for our 
constraints. 
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FIG. 13: The asymmetric curve corresponds to a local maximum for all points on the segment from A (where it crosses the 
line ^ — rj + 4/3/3^) to B (where it crosses ^ — — 4//3^ and becomes unphysical. It is the global maximum only in the region 
below the diagonal, which it crosses at D. If /3 < 2, _D lies above B so the asymmetric configuration is never a global maximum. 
If /9 > 2 then the asymmetric configuration is the global maximum for every point on the segment between D and B. 



join. The bottom of the lens meets the diagonal at the (3 = 2 line. 
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